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(This  Abstract  is  Unclassified) 


ABSTRACT 


Analytical  studies  to  determine  the  characteristics  of  textile  para¬ 
chute  canopies  in  supersonic  flow  and  associated  airflow  phenomena  are 
described  to  the  point  of  completion  in  this  program. 

f”  '  ” 

Data  derived  from  the  wnnd  tunnel  test  progran  ,,  which  v/ere  conduc¬ 
ted  using  six  inch  and  eight  inch  parachute  models,  are  cJijJcussed.  Test 
Mach  numbers  covered  the  range  from"  with  approxi¬ 

mately  200  test  conditions  tabulated.  The  over-all  test  program  discussion 
describes  investigations  leading  to  the  evolution  of  canopy  shapes  which 
performed  satisfactorily  through  Mach  4.  65. 

Force  data  and  high  speed  Schlieren  photographs  of  numerous  test 
conditions  are  included  to  illustrate  canopy  performance  characteristics, 

-  Total  canopy  porosities  of  8  to  18  percent  with  high  roof  porosities 
and  low  inlet  region  porosity  distributions  were  studied  most  intensively. 


Canopy  shapes  consisting  of  bi-conic,  conic,  modified  guide  surface, 
and  conventional  ribbon  types  were  investigated. 
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SECTION  I 


STATUS  OF  PARACHUTE  DEVELOPMENT  AT  SUPERSONIC  SPEEDS 
A.  Introduction 

At  the  time  of  the  initiation  of  this  program,  June  30,  1961,  a  consider¬ 
able  amount  of  parachute  testing  had  been  performed  at  transonic  and  super¬ 
sonic  speeds.  These  test  programs  were  conducted  using  high  speed  sleds, 
wind  tunnels,  and  free  flight  vehicles.  A  large  amount  of  data  was  generated 
as  a  result  of  these  programs  -  such  as  force  measurements,  high  speed  pho¬ 
tographs,  and  Schlicren  films.  Numerous  conventional  canopy  configurations 
were  employed  in  these  tests  with  ribbon-type  constructions  predominating. 

In  this  part  of  the  report,  a  brief  background  of  parachute  testing  and  develop¬ 
ment  will  be  presented.  As  a  result  of  analysis  of  the  data  derived  from  the 
early  tests,  the  basic  problems  associated  with  the  application  of  conventional 
parachute  configurations  at  supersonic  speeds  were  determined  and  will  be 
discussed.  Section  I  will  be  devoted  to  this.  Subsequent  sections  of  this  report 
will  describe  the  objectives,  procedures  and  results  of  this  program. 

This  program  was  initiated  in  6rder  to  explore  the  applicability  of 
flexible,  self-inflating  parachute -like  configurations  in  the  Mach  number 
range  from  M  =  1.  5  to  M  =  4.  0.  The  program  was  basically  divided  into  two 
separate,  concurrent,  yet  interdependent  efforts.  These  consisted  of:  (1) 
a  theoretical  study,  and  (2)  a  wind  tunnel  test  program.  The  theoretical 
study  had  as  its  objective,  the  development  of  a  mechanized  rigorous  theo¬ 
retical  method  of  prediction  of  the  flow  field  about  a  trailing  decelerator 
leading  to  the  determination  of  the  surface  properties  of  the  decelerator. 

The  development  of  such  a  computer  program  serves  to  define  the  environ¬ 
ment  surrounding  a  parachute-like  configuration  sc  that  a  shape  which  will 
perform  satisfactorily  in  stability  and  inflation  characteristics  at  supersonic 
speeds  may  be  predicted.  The  effects  of  a  non-uniform  forebody  wake  was 
also  considered  in  the  formulation.  The  wind  tunnel  program  was  initially 
conceived  to  be  a  testing  program  which  verified  the  results  of  the  theoreti¬ 
cal  prediction.  Due,  however,  to  formulation  and  computer  programming  diffi¬ 
culties  encountered  in  the  analytical  study,  this  program  was  unable  to  pro¬ 
vide  such  predictions  consistent  with  available  Government  facility  wind  tunnel 
testing  dates.  The  prediction  of  feasible  canopy  configurations  was  performed 
by  non-rigorous  procedures.  These  involved  the  study ofSchlierenphotographs 
of  early  tests,  coupled  with  the  application  of  shocktables  and  employing  an 
empirical  formulation  of  wake  flow  properties.  This  procedure  led  to  the  pre¬ 
diction  of  initial  modifications  to  conventional  parachute  configurations  which 
were  indicated  to  be  required  to  improve  their  performance  characteristics. 
Successive  testing  programs,  in  conjunction  with  this  procedure,  resulted 
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in  the  gradual  improverneni  of  canopy  performance  and  led  to  the  catabliah- 
ment  of  the  apparently  significant  phenomena  which  defines  the  criteria  for 
the  determination  of  satisfactory  parachute  performance.  Unconventional 
parachute  canopy  designs  which  successfully  endured  tests  at  Mach  numbers 
through  M  =4.65  were  developed  in  this  program.  The  success  of  these 
canopies  in  the  final  (Phase  II)  testing  program  resulted  in  security  classifi¬ 
cation  of  the  designs  of  these  configurations. 

B.  Basic  Problems  Associated  with  Supersonic  Operation  of  Conventional 
Parachute  Configurations 

As  discussed  in  Subsection  A  above,  numerous  test  programs  have  been 
accomplished  during  the  investigation  of  the  performance  of  conventional  par¬ 
achutes  at  supersonic  speeds.  However,  such  configurations  have  been  found 
to  perform  poorly  in  the  supersonic  speed  range,  particularly  for  Mach  num¬ 
bers  above  1.  8  to  2,  0.  Specifications  for  a  particular  decelerator  configura¬ 
tion  require,  in  general,  that  the  ratio  of  drag  effectiveness  to  system  weight 
be  a  maximum,  and  further  that  the  voliunetric  requirements  of  the  system 
be  minimized.  When  self-inflating  parachute  configurations  properly  operate, 
the  system  drag  to  weight  ratio  is  lower  than  other  decelerator  types  such  as 
inflatable  devices.  For  example,  an  inflatable  drag  device,  whether  spheri¬ 
cal  or  conical  in  shape  represents  a  significantly  greater  total  material  sur¬ 
face  area  (hence  weight)  for  a  given  drag  area.  A  drag  coefficient  of  approxi¬ 
mately  0,  2  based  upon  total  cloth  area  is  typical  for  an  inflatable  device  (not 
including  inflation  aids,  etc.,)  whereas  parachute  configurations  have  achieved 
drag  coefficients  approaching  1,  0  based  upon  cloth  area  in  the  Mach  number 
range  for  which  efficient  operation  of  the  conventional  parachute  is  realized. 
However,  for  Mach  numbers  above  1.  8  to  2.  0,  the  conventional  parachute 
drag  efficiency  has  been  found  genercdly  to  decrease  rapidly  with  increasing 
Mach  number.  In  fact,  the  drag  effectiveness  of  conventional  parachutes  has 
been  found  to  be  in  the  order  of  1  /6  of  the  expected  value  at  a  Mach  number  of 
3.  5.  This  would  infer  a  cross-over  in  drag  efficiency  of  an  inflatable  device, 
for  illustration,  as  compared  to  a  conventional  parachute  (i,  e, ,  previously 
employed  configurations)  for  Mach  numbers  in  excess  of  this  value. 

The  following  paragraphs  of  this  section  will  describe  the  characteris¬ 
tics  of  conventional  parachute  configurations  when  operating  in  the  supersonic 
range  and  contribute  to  their  reducing  efficiency  in  an  increasing  Mach  num¬ 
ber  environment. 

1.  Inflation  Characteristics 

Analysis  of  schlieren  movies  of  conventional  parachute  configura¬ 
tions  has  revealed  a  typical  behavior  in  the  supersonic  range  above 
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Mach  1.8  -  2,0.  This  behavior  is  characterized  by  a  rather  violent 
breathing  tendency,  or  alternate  inflation  and  deflalion  cycle,  which 
occurs  at  very  high  frequencies  {the  order  of  50  to  100  cycles  per 
second).  This  phenomena  is  illustrated  in  Figure  I-l  for  a  4,  0  foot 
diameter  flat  circular  ribbon  parachute  with  a  19  percent  porosity. 

This  film  sequence  was  taFen  during  a  wind  tunnel  test  of  this  configura¬ 
tion  at  a  Mach  number  of  3.  5,  during  the  conduct  of  a  series  of  tests 
in  the  NASA,  Cl  e '/eland,  Ohio,  10  feet  by  10  feet  section  by  these 
Laboratories  under  Contract  AF  33(6l6)-3346.  During  this  program, 
both  ribbon  and  guide  surface  type  parachutes  were  tested  in  the  Mach 
number  range  from  M  =  2.  0  to  3.  5.  A  wide  number  of  variations  of 
both  types  of  parachutes  was  employed.  For  example,  the  porosity  of 
the  ribbon-type  parachutes  was  varied  from  5  to  30  percent.  Porosity 
distributions  were  also  considered  in  that  low  skirt  porosity  with  large 
vent  porosity  as  well  as  low  vent  porosity  and  high  skirt  porosities 
were  employed  in  the  various  configurations.  In  addition,  the  tests  were 
conducted  both  with  and  without  a  fore  body.  In  all  instances,  parachute 
performances  from  the  standpoint  of  inflation  stability  were  as  mention¬ 
ed  above  wdth  high  frequency  breathing  evident. 

In  addition  to  the  ribbon  parachutes,  various  t'ypes  of  Guide  Sur¬ 
face  Ribless  canopies  were  also  tested.  The  extremes  in  breathing  are 
not  as  pronounced  as  for  the  ribbon;  however,  overinflation  of  the  roof 
and  indentation  of  the  guide  surfaces  resulted  in  nearly  immediate  fail¬ 
ure  of  the  canopy  at  the  intersection  of  the  roof  and  guide  surface  panels. 

In  order  to  evaluate  whether  this  phenomena  was  associated  with 
the  flexible  materials  of  which  these  configurations  are  constructed  or 
to  an  unstable  flow  field  aJ.ead  of  the  canopy,  a  further  test  program 
was  conducted  employing  rigid  canopies  constructed  of  steel  and  simu¬ 
lating  the  inflated  shape  of  a  typical  ribbon  canopy.  These  tests  were 
conducted  with  and  without  simulated  lines,  various  numbers  of  lines, 
and  lengths  of  lines.  Based  upon  the  results  of  these  tests,  the  exist¬ 
ence  of  unsteady  flow  conditions,  and  Mach  number  effects  was  con¬ 
firmed.  A  discussion  of  these  tests,  and  the  preceding  fabric  model 
tests  along  with  an  analysis  of  the  various  associated  flow  regimes  is 
given  in  Reference  1.  As  a  result  of  this  analysis  it  appeared  that  the 
nature  of  the  flow  regime  and  shock  wave  structure  ahead  of  the  canopy 
is  significantly  dependent  upon  the  forebody  wake  characteristics,  and 
that  parachute  configurations  are  greatly  affected  by  the  interaction  of 
the  canopy  shock  wave  with  velocity  gradients  existing  in  the  wake  of  a 
forebody.  Several  typical  flow  types,  which  are  associated  with  this 
interaction,  are  discussed  in  Reference  1. 
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Figure  I-l.  Characteristics  of  Two  Parachutes  at  M  =  3.  5 
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In  later  tests,  a  canopy  confi{>nration  which  was  tested  in  the 
Cook  Technological  Center  wind  tunnel  facility  (under  Contract 
AF  3i(6l6)-5507),  demonstrated  that  certain  configurational  charac¬ 
teristics  may  have  a  strong  effect  in  suppressing  the  unfavorable 
interaction  effects  which  appear  to  be  associated  with  a  parachute  con¬ 
figuration  immersed  in  a  non-uniform  flow  field  provided  by  an  up¬ 
stream  forebody.  The  canopy  of  this  configuration  consisted  of  a  15 
(half  angle)  conir.iJ  frustrum  fabricated  of  a  non-porous  material  with 
a  vent  which  pro\'ided  an  exit  to  inlet  area  ratio  of  0.  Z9.  A  variation 
of  this  configuration  employing  a  45°  conical  frustrum  was  also  tested. 
Very  high  drag  coefficients  (0.  •15  and  0.  96,  respectively)  based  upon 
cloth  area,  were  realized  with  these  models  at  M  =  2.  Extremely  high 
inflation  characteristics  were  also  evident  from  the  schlieren  films 
(Figure  1-2)  and  very  little  canopy  breathing  was  evident.  How-ever, 
these  canopies  exhibited  an  instability  about  the  point  of  suspension  by 
oscillation  angles  up  to  nearly  20°.  In  spite  of  the  high  oscillation 
angles,  the  canopies  remained  well  inflated  although  subjected  to 
severe  disturbances  which  were  associated  with  the  high  angular  ex¬ 
cursions,  The  instability  of  these  configurations  was  attributed  to  the 
ring  airfoil-type  lift  inherent  in  this  canopy  shape,  A  stabilizing  sur¬ 
face  extending  forw-ard  of  the  maximum  diameter  had  been  added  to  the 
45*^  conical  canopy.  However,  the  high  inflation  tendencies  of  this 
canopy  resulted  in  material  stretch  of  this  surface  to  such  an  extent 
that  it  was  ineffective  in  contributing  st.abilization,  It  was  this  basic 
configuration,  along  with  the  indicated  necessary  modifications  which 
served  to  establish  some  of  the  configurational  details  of  the  shapes 
which  were  selected  for  the  first  testing  series  in  this  program, 

2,  Stability  Characteristics 

Analysis  of  the  data  from  early  tests  indicate  that  the  stabilityof 
a  parachute  canopy  about  the  point  of  suspension  varies  directly  with 
porosity.  Tests  in  the  present  program  indicate  this  trend  also.  How¬ 
ever,  it  is  apparent  that  the  configurational  details  of  the  canopy  have 
a  significant  effect  also.  An  illustration  of  this,  of  course,  is  the 
guide  surface  type  canopy  which  is  constructed  generally  with  a  low 
porosity  material  yet,  in  the  subsonic  and  transonic  ranges  of  Mach 
numbers,  exhibits  extreme  stability.  The  contribution  of  the  forward 
guide  surface  is  considered  significant  in  the  attainment  of  this  stabil¬ 
ity.  The  flat  roof  of  this  configuration  contributes  negative  lift  when 
sufficient  porosity  allows  flow  through  the  canopy.  The  forward  guide 
surfaces,  which  are  inflated  internally  with  near- stagnation  pressures, 
serve  as  a  stabilizing  cone-frustrum  and  contribute  the  necessary  com¬ 
pensating  positive  lift  for  stability.  The  basic  guide  surface  configura¬ 
tion,  along  with  modifications,  served  also  to  provide  some  of  the 
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shapes  investigated  in  the  first  testing  series  of  this  program  These 
configurations  will  be  discussed  in  Section  II, 

3.  Shock  Wave-Boundary  Layer  Interaction 

Analysis  ofSchlieren  photographs  of  conventional  parachute  con¬ 
figurations  at  supersonic  speeds  has  indicated  significant  shock  wave- 
suspension  line  boundary  layer  interactions.  It  has  appeared  from 
observation  of  parachute  test  data,  that  for  those  configurations  which 
exhibited  low  or  marginal  inflation  characteristics  that  the  effect  of 
the  shock  wave-boundary  layer  interaction  has  a  significant  effect  on 
the  inflation  stability  of  the  canopy,  and  has  appeared  to  be  somewhat 
associated  with  cyclic  canopy  collapse.  In  other  instances  (Figure  1-2) 
it  has  appeared  that  this  phenomenon  does  not  have  a  significant  effect 
on  the  parachute  performance.  This  latter  example,  of  course,  applied 
to  a  canopy  with  strong  inflation  properties  as  previously  described. 

This  phenomenon  has  been  most  evident  in  those  instances  where 
the  suspension  lines  become  normal,  or  near  normal  to  the  canopy 
shock  wave.  This  condition  is  readily  attained  in  the  case  of  a  canopy 
which  is  unstable  about  the  point  of  suspension,  allowing  excursion 
angles  of  the  canopy  such  that  the  lines  approach  normalcy  to  the 
shock  wave,  such  as  illustrated  in  Figure  1-2,  This  is  possible,  since 
examination  of  Schlieren  films  appears  to  indicate  that  the  canopy 
shock  wave  translates  normal  to  the  wake  axis,  an  amount  approximate¬ 
ly  equed  to  the  width  of  the  fore  body  wake.  Very  small  rotation  of  the 
shock  pattern  about  its  vertex  was  evident.  For  those  canopy  types 
which  sustain  a  strong  shock,  very  small  oscillation  angles  are  needed 
to  attain  a  near  normal  angle  between  the  suspension  lines  and  the 
shock  wave. 

Since  separation  of  a  boundary  layer  by  shock  wave  interaction 
is  associated  with  a  critical  pressure  ratio  (as  a  function  of  Mach 
number)  across  the  shock  in  the  case  of  a  shock  wave  incident  to  a 
flat  plate  (Reference  2),  it  is  expected  that  this  will  be  true  also  for  a 
shock  wave  intersection  with  a  three  dimensional  body.  The  Basic 
Test  program  was  organized  to  establish  the  effects  of  shock  wave 
interaction  with  a  boundary  layer  on  a  cylinder  as  a  function  of  shock 
wave  angle,  shock  wave-cylinder  intersection  angle,  Reynolds  num¬ 
ber  and  Mach  number.  The  Basic  Test  program  will  be  discussed  in 
Section  II. 

C,  Indications  of  Desirable  Canopy  Geometry  Based  Upon  Early  Tests 
Reference  was  made  in  Subsection  B  of  Section  I  to  a  conical  canopy 
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configuration,  illustrated  in  Figure  I-Z  which  had  exhibited  strong  inflation 
tendencies  as  well  as  inflation  stability  at  Mach  2.  In  the  design  of  this 
canopy,  three  features  differed  extremely  from  conventional  configurations 
previously  applied  in  the  supersonic  speed  range.  These  are  (1)  a  much 
lower  exit  to  inlet  area  ratio  had  been  employed,  (2)  a  non-porous  material 
was  used  in  the  vicinHy  of  the  inlet,  or  skirt  region  of  the  canopy,  and  (3) 
a  much  lower  total  porosity  was  employed. 

The  total  porosity  of  a  parachute  is  defined  as  the  open  area  divided  by 
the  total  cloth  area.  Thus,  the  open  area  consists  of  the  space  between  rib¬ 
bons,  as  in  a  ribbon  canopy  construction  and  open  spaces  in  a  woven  fabric 
which  allow  the  passage  of  air.  In  the  case  of  a  canopy  such  as  described 
above  which,  except  for  the  open  vent,  is  constructed  of  non-porous  material, 
the  over-all  porosity  was  considered  as  the  area  of  the  vent  divided  by  the 
total  area  (which  includes  the  vent).  This  definition  will  be  at  variance  to 
some  degree  with  porosity  which  is  expressed  as  a  function  of  permeability. 
The  latter  property  of  cloth  is  determined  by  measuring  the  rate  of  flow  of 
air  under  specified  conditions  through  a  fabric.  The  normally  small  spaces 
between  fibers  will  prescribe  somewhat  different  flow  conditions  than  flow 
through  a  large  vent.  However,  in  all  calculations  performed  in  this  pro¬ 
gram  to  determine  the  open  area,  a  strictly  geometric  interpretation  has 
been  made.  This  "open"  area  is  used  to  determine  both  totad  porosity  and 
"exit"  area. 

Accordingly,  the  exit  to  inlet  area  ratio  of  the  conical  canopy  discussed 
above  is  0.  29,  and  the  total  geometric  porosity  is  9.  5  percent.  In  contrast 
to  this,  a  conventional  ribbon  parachute  is  constructed  using  a  geometric 
porosity  of  1  5  to  28  percent.  If  the  canopy  were  constructed  in  a  hemispheri¬ 
cal  shape,  the  exit  to  inlet  area  ratio  would  be  0.  30  to  0,  56,  assuming  full 
inflation.  In  such  configurations,  the  porosity  at  the  skirt  is  normally 
greater  than  that  near  the  vent,  making  this  area  less  able  to  sustain  the 
near-stagnation  pressures  expected  in  the  canopy  inlet  region.  Strong  infla¬ 
tion  characteristics  in  the  skirt,  or  inlet  regions  appear  to  be  essential  to 
provide  an  outward  force  at  least  equal  to  the  inward  components  of  the  sus¬ 
pension  line  loads. 

As  discussed  previously,  the  stability  of  a  given  parachute  configura¬ 
tion  about  the  point  of  suspension  appears  to  vary  inversely  with  porosity. 

In  the  case  of  the  Guide  Surface  configuration,  however,  the  total  porosity 
is  low  and  the  stability  is  good,  although  structural  damage  resulted  in  tests 
of  this  configuration  at  supersonic  speeds.  The  stabilizing  effect  of  the  for¬ 
ward  guide  su^'face  is  apparent,  although  the  (conventionally  employed)  angle 
of  this  surface  appeared  to  be  too  great. 
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The  apparent  requirements  for  low  over -all  porosity,  coupled  with 
the  use  of  stabilizing  surfaces  appeared  essential.  A  study  of  the  effects  of 
(1)  porosity  distributions  in  producing  good  inflation,  (2)  stabilizing  surface 
angles  and  lengths  for  proper  inflation  and  stability,  and  (3)  the  over-all 
geometric  shape  of  the  canopy  for  optimum  drag  to  weight  ratio  appeared  to 
prescribe  the  direction  of  effort  in  the  Phase  1  parachute  test  program. 

D.  Effects  of  Forebody  Wake  on  Supersonic  Parachute  Operation 

The  importance  of  the  effects  of  the  wake  of  a  forebody  on  a  trailing 
parachute  configuration  have  been  recognized.  In  Reference  1,  calculations 
were  performed  to  determine  the  effects  of  the  non-uniform  flow  field  pro¬ 
duced  by  the  wake  on  a  shock  wave  ahead  of  a  blunt  trailing  decelerator. 
Significant  shock  diffraction  was  verified  by  employing  empirical  wake 
theory  for  the  determination  of  the  wake  flow  ahead  of  the  canopy  and  shock 
wave.  Diffraction  of  the  typical  near-normal  shock  wave,  which  would 
exist  under  free  stream  conditions  ahead  of  a  blunt  body  such  as  an  inflated 
parachute,  into  a  conical  type  shock  was  indicated  by  these  calculations. 

This  has  been  confirmed  by  numerous  Schlieren  photographs  of  parachutes 
operating  supersonically  in  a  forebody  wake.  Furthermore,  the  dependence 
upon  the  wake  of  a  forebody  for  satisfactory  parachute  operation  at  super¬ 
sonic  speeds  has  also  been  indicated  by  past  tests  and  applications. 

Examination  of  a  Phase  la  model  configuration  under  test  will  be  used 
as  an  illustration  to  establish  the  flow  field  ahead  of  a  canopy  and  associated 
shock  wave.  To  predict  the  flow  field  in  the  wake,  an  empirical  wake  theory 
will  be  employed.  Drag  predictions  provided  by  this  procedure  on  trailing 
decelerators  have  been  in  good  agreement  with  test  data. 

The  velocity,  V^,  across  the  wake  as  given  in  Reference  1: 


where:  ' 

y  =  the  distance  from  the  wake  centerline 
b  =  the  wake  half  width 
Vqo  =  the  free  stream  velocity 

^Wi 

_ :  =  the  wake  centerline  velocity  decrement  ratio  which  is 
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determined  by  the  expression 


^Wi  /^oo 


C 


(x/d  -  c 


where  the  constant  C  is  based  upon  test  data.  (0.  7  used  in  this  report) 

Figure  1-3  is  a  trace  from  a  Schlieren  photograph  of  a  Phase  la 
configuration(discussed  in  Section  U),  which  was  located  at  an  x/d  =  10,  at  a 
test  Mach  number  of  2.  13.  Various  wake  Mach  numbers  are  indicated  In  the 
wake  ahead  of  the  shock  wave  and  flow  deviations  through  the  shock  for  these 
Mach  numbers  and  local  shock  wave  angles  are  also  shown.  A  wake  width 
equal  to  tne  forebody  diameter  is  assumed.  It  may  be  seen  from  Figure  1-3 
that  a  subsonic  outflow,  resulting  from  the  deflection  through  the  shock  has 
a  very  favorable  effect  in  maintaining  the  canopy  fully  inflated.  Canopy 
inlet  angles,  below  the  horizontal,  approaching  25  degrees  will  provide  a 
positive  "angle  of  attack"  of  the  canopy  inlet  to  the  local  flow  direction.  It  is 
apparent  from  Figure  1-3  that  the  flow  deviation  angle  based  upon  the  free 
stream  Mach  number  is  sufficient  to  spread  the  wake  to  near  the  skirt  of  the 
canopy.  This  is  further  indicated  to  be  required,  since  the  pressure  ratio 
across  the  normal  portion  of  the  shock  on  the  centerline  is  greater  (for  the 
wake  centerline  Mach  number)  than  the  pressure  ratio  across  the  oblique 
shock  at  the  free  stream  Mach  number.  Divergence  of  the  wake  flow  is 
further  encouraged  by  the  high  (near  stagnation)  pressures  existing  in  the 
canopy.  Since  wake  effects  of  the  parachute  confluence  point  are  not  consid¬ 
ered,  it  is  expected  that  the  wake  Mach  numbers  near  the  centerline  are 
lower  than  predicted  giving  a  lower  than  predicted  pressure  ratio  across  the 
shock  at  the  centerline  such  that  a  near  constant  pressure  distribution  is 
realized  throughout  the  diverged  wake  ahead  of  the  canopy,  satisfying  pres¬ 
sure  boundary  requirements. 

In  the  event  that  the  parachute  diameter  ratio  to  the  wake  width  were 
significantly  greater,  and  assuming  a  similar  shock  geometry,  the  flow 
deviation  through  the  shock  wave  would  not  be  adequate  to  spread  the  wake 
(if  the  shock  standoff  distance  remained  the  same)  to  the  skirt  of  the  canopy, 
and  the  flow  aft  of  the  shock  at  greater  distances  from  the  wake  centerline 
would  be  supersonic.  A  lip  shock  would  then  be  produced  which  may,  as 
past  tests  indicate,  further  spread  the  wake  upon  intersecting  the  boundary 
layer  until  the  entire  canopy  is  in  subsonic  flow.  This  would  possibly  lead 
to  a  fluctuating  flow  condition,  and  reduction  in  parachute  performance.  In 
consideration  of  this,  and  other  associated  phenomena,  it  is  expected  that 
a  practical  parachute  to  forebody  diameter  ratio  limit  may  possibly  exist 
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for  efficient  supersonic  parachute  operation  at  a  prescribed  supersonic 
Mach  number, 

E.  Purpose  and  Scope  of  This  Program 

As  pointed  out  previously,  a  number  of  appj  ied  research  programs  have 
been  conducted  to  study  and  to  extend  the  operational  capabilities  of  textile 
parachute  canopies  into  the  supersonic  speed  regime.  Results  of  these  pro¬ 
grams  revealed  problem  areas  which  required  further  study  in  order  to 
evolve  satisfactory  concepts  of  self- inQaling  aerodynamic  decelerators  for 
supersonic  speed  applications. 

This  program  was  planned  to  consist  of  an  exclusive  analysis,  wind 
tunnel  testing,  data  documentation,  and  reporting  of  phenomena  associated 
with  the  operation  of  flexible  type  self-inflating  aerodynamic  decelerators 
operating  within  a  Mach  number  of  range  of  1.  5  to  4.  0. 

The  analytical  study  was  to  consist  of  an  application  of  fluid  flow 
equations  to  configurations  consisting  of  parachute-like  shapes  with  flow 
properties  surrounding  these  configurations  described  by  means  of  the  inte¬ 
gration  of  established  fluid  flow  equations.  Mechanization  of  all  computa¬ 
tions  was  to  be  maximized,  and  correlation  of  derived  fluid  flow  properties 
with  wind  tunnel  observations  was  to  serve  (1)  to  establish  the  adequecy  of 
developed  techniques,  or  (2)  to  serve  as  a  basis  for  modification  of  tech¬ 
niques,  so  as  to  produce  the  required  degree  of  agreement. 

The  ultimate  goal  of  the  analytical  approach  was  to  evolve  an. aerody¬ 
namic  theory  which  would  serve  to  both  predict  and  explain  the  test  results, 
and  lead  to  a  means  of  prediction  of  parachute  canopy  geometry  for  satis¬ 
factory  performance  at  supersonic  speeds. 

An  exploratory  wind  tunnel  program  was- designed  to  be  conducted  in 
order  to  supplement  and  verify  the  analytical  study. 

In  subsequent  sections  of  this  report,  the  analytical  procedures  and 
exploratory  test  programs  will  be  described. 
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SECTION  U 


EXPERIMENTAL  PROGRAM 


Three  test  phases  were  completed  in  this  program.  The  Phase  I  test 
program  was  conducted  In  the  Cook  Technological  Center  wind  tunnel  and  the 
AEDC-VKF,  Tunnel  A,  Tullahoma,  Tennessee.  The  Interim  test  program 
(Phase  la)  was  conducted  In  the  Cook  Technological  Center  wind  tunnel  facility, 
and  the  third  program  (Phase  U)  was  conducted  in  the  Langley  Research  Center 
Large  Supersonic  Tunnel,  Langley  Field,  Virginia.  In  the  three  phase  test 
program,  tests  were  conducted  over  a  Math  number  range  of  from  1.  5  to  4.65. 

In  order  to  make  maximum  use  of  available  testing  time  in  Government 
wind  tunnel  facilities,  the  first  testing  program  was  required  to  be  conducted 
earlier  in  the  over-all  program  than  initially  planned.  As  a  result,  this  pro¬ 
gram  was  required  to  be  of  a  development  nature,  and  model  designs  were 
accordingly  required  to  be  determined  by  general  trends  in  geometric  shapes 
which  were  indicated  to  be  favorable  in  earlier  testing  programs,  as  dis¬ 
cussed  in  Section  I,  Parts  C  and  D, 

A.  Phase  I  Teat  Program 

1.  Phase  I  Test  Facilities 

The  first  test  phase  of  this  program  was  conducted  at  the  Cook 
Technological  wind  tunnel  facility  and  the  Arnold  Engineering  Develop¬ 
ment  Center,  Tullahoma,  Tennessee.  The  AEDC  tests  were  conducted 
in  the  40  inch  x  40  inch  tunnel  A  test  section  which  is  capable  of  a  Mach 
number  range  of  1.  5  to  5.  A  description  of  these  test  facilities  is  given 
in  Appendices  I  and  11. 

2.  Basic  Test  Program 
a.  General 

As  a  result  of  the  study  of  various  Schlieren  photographs 
of  past  tests,  shock  wave-suspension  line  boundary  layer  inter¬ 
action  appeared  to  be  evident  as  discussed  in  Section  I.  A  con¬ 
siderable  amount  of  data  appears  to  be  available  in  the  literature 
on  two-dimensional  shock  wave-boundary  layer  interaction. 
Pressure  ratios  associated  with  separation  of  the  boundary  layer 
on  a  flat  plate  are  fairly  well  established  (Ref.  2),  and  a  certain 
degree  of  correlation  of  this  data  is  associated  with  parachute 
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data  (Ref.  1). 


In  order  to  eatabliali  correlation  of  publiahed  two-dimension¬ 
al  data  with  three-dimensional  characteristics,  the  basic  test  pro¬ 
gram  of  this  contract  was  conceived;  it  was  decided  to  employ 
rigid  models.  Since,  as  discussed  above,  critical  interaction  con¬ 
ditions  appear  to  exist  when  the  shock  wave  and  suspension  line 
are  approximately  normal  to  each  other,  appropriate  test  items 
appeared  to  consist  of  a  shock  generator  to  provide  a  nearly  nor¬ 
mal  shock,  and  a  rod  to  simulate  a  line.  An  adjustment  mechan¬ 
ism  would  provide  for  a  variation  of  angle  between  the  shock  wave 
and  the  rod,  by  variation  of  the  rod  angle  of  attack.  The  effect 
of  Reynolds  number  would  be  determined  by  an  axial  position  ad¬ 
justment  of  the  length  of  the  rod  ahead  of  the  shock  wave  or  by 
variation  of  tunnel  total  pressure,  /.r.  apparently  appropriate  con¬ 
figurational  design  was  conceived  for  preliminary  tests  in  the  CTC 
wind  tunnel  facility.  This  is  shouTi  in  Figure  II-l,  The  shock 
generator  employed  in  this  test  configuration  consisted  of  a  solid 
choked  annular  ring  allow'ing  flow  through  the  center  of  the  ring, 
and  also  aUowing  freedom  of  motion  of  the  rod  inserted  through  it 
and  penetrating  the  shock  generated  by  the  choked  annulus. 

b.  Cook  Technological  Center  Basic  Test  Program 

These  tests  were  conducted  by  varying  the  rod  angle  every 
2  degrees  from  0°  to  22°,  and  at  each  angular  position  a 
Schlieren  photograph  was  taken.  Table  II-l  describes  the  test 
schedule.  This  procedure  was  repeated  providing  three  different 
rod  lengths  ahead  of  the  shock  wave.  Analysis  of  the  Schlieren 
photographs  indicated  that  precise  measurements,  the  nature  of 
which  were  desired  in  these  tests,  could  not  be  performed  with 
the  test  configuration  employed.  An  inherently  unstable  flow 
condition  had  been  encountered  since  upon  establishment  of  a 
normal  shock  in  front  of  the  annulus,  and  resultant  separation  of 
the  boundary  layer  (Figure  II-2)  the  separated  flow  is  not  stagna¬ 
ted  on  the  blunt  nose  (unlike  separation  ahead  of  a  blunt  body  by 
means  of  a  protruding  spike).  Once  separation  is  induced,  the 
separated  flow  must  feed  through  the  annulus,  thereby  establish¬ 
ing  the  need  for  a  new  normal  shock  on  the  axis  to  be  established 
to  maintain  mass  flow  requirements  through  the  annulus.  This 
apparently  occurs  at  a  high  frequency,  and  the  camera  speeds 
employed  were  not  adequate  to  stop  this  high  frequency  phenom¬ 
enon.  The  tests  indicated  that  in  order  to  isolate  the  phenomena 
of  interest,  the  rod-shock  wave  intersection  must  be  outside  of 
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Figure  II-2  Typical  Schlieren  Photograph,  Phase  I  Basic  Test,  CTC  Wind 
Tunnel.  Film  Speed;  E4  fps,  M  =  2.  1,  a  =  8° 
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the  ring,  so  as  to  be  unaffected  by  the  downstream  flowthrough 
the  annulus.  In  anticipation  of  a  more  thorough  and  detailed 
study,  the  test  tnethod  was  modified  for  the  basic  test  phase  in 
the  AEDC  Tunnel  A.  This  test  progr;un  and  associated  equip¬ 
ment  are  discussed  below, 

c.  AEDC  Basic  Test  Program 

In  Figure  II- 3  is  shown  the  test  installation  in  the  AEDC 
Tunnel  A,  where  the  investigation  of  shock  wave-boundary  layer 
interaction  studies  were  performed  as  part  of  the  Phase  I  test 
program. 

A  remotely  controlled  actuator  which  is  capable  of  limited 
transverse  and  angular  movement,  to  which  the  rod  is  attached, 
coupled  with  the  AEDC  Tunnel  A  positioner  to  which  the  shock 
generator  is  attached,  provides  a  system  v/hich  is  capable  of 
producing  relative  variation  in  two  linear  dimensions  and  one 
angular  (pitch)  movement  for  the  rod  mounted  to  the  actuator, 
and  the  shock  generator  w’hich  is  mounted  to  the  axial  positioner. 

This  investigation  w-as  conducted  over  a  Mach  number  range 
of  from  M  =  1.  5  to  M  =  3,  with  variation  of  Reynolds  number  in  < 

some  cases.  The  test  program  is  described  in  Table  II-Z,  In 
Figure  LI-4  is  shown  the  shadowgraph  pictures  of  the  shock  wave 
produced  by  the  shock  generator.  These  photographs  were  made 
by  focusing  the  Schlieren  system  on  the  rear  wind  tunnel  window. 

Traces  of  these  shock  wave  configurations  were  made  during 
wind  tunnel  test  runs  which  were  performed  during  Schlieren 
system  checkout.  Some  considerable  degree  of  difficulty  with  the 
Tunnel  A  Schlieren  system  was  encountered  and  the  shock  wave 
sketches  above  were  made  in  order  to  save  tunnel  running  time 
in  accomplishing  basic  tests.  The  test  rod  (Figure  II-3)  is 
equippped  with  two  sets  of  four  holes  each.  The  sets  are  located 
an  inch  apart,  with  holes  in  each  set  located  at  90°  intervals 
around  the  rod  and  in  a  plane  perpendicular  to  the  rod  axis.  From 
these  holes,  pressure  lines  were  run  through  the  rod  and  axial 
positioner  hub  to  the  outside  of  the  tunnel.  Static  pressure 
measurements  were  thereby  made  possible  at  desired  rod  posi¬ 
tions.  The  rod  was  positioned  such  that  the  shock  wave  location 
would  be  intermediate  to  the  two  sets  of  static  orifices.  Since 
the  shock  wave  position  could  not  be  adequately  determined  » 

theoretically,  the  shadowgraphs  described  above  were  used. 

Accordingly,  the  coordinates  of  three  points  on  the  shock  wave 
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Figure  II-3.  Phase  I  Basic  Test  Installation 
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igiiT6  XI  4.  Typical  Shock  Con£igu.ra.tioris  X^ha.S0  1  S3.sic  Tests 


(at  whicli  tlu'ff  local  shock  wave  angles  were  incasurcd)  were 
deterinined  relative  to  the  control  counts  oti  the  actuator  and  on 
the  tunnel  axial  positioner.  On  the  basis  of  these  initial  posi¬ 
tions,  actuation  niechanism  positions  were  then  determined  so 
as  to  provide  an^idar  variatiori  of  the  rod  relative  to  the  shock 
wave,  while  niaintainin^  the  shock  wave  intermediate  to  the  two 
sets  of  static  pressure  taps  (providine  a  constant  rod  length 
ahead  of  the  shock  v..ive).  It  liiay  be  added  liere  that  since  the 
rod  actuation  mechanism  was  not  capable  of  axial  movenicnt, 
movement  of  the  rod  in  an  axial  direction  relative  to  the  shock 
wave  was  accomplisht;d  by  axial  movement  of  the  shock  genera¬ 
tor  (provided  by  the  tunnel  axial  positioner  on  which  it  was 
mounted),  resulting  in  axial  translation  of  the  shock  wave. 

The  tests  as  outlined  in  Table  il  -2  were  conducted  in  sub¬ 
sequent  tunnel  A  runs.  The  correlation  of  rod  angle,  inception 
of  boundary  layer  separation  on  the  rod  and  associated  pressure 
measurements  before  and  after  separation  were  dependent  upon 
the  inspection  of  the  corresponding  still  Schlieren  films  and 
Schlieren  movies.  Due  to  the  difficulty  with  the  tunnel  Schlieren 
system,  a  negligible  number  of  Schlieren  pictures  were  provided. 
Consequently  the  tests  contributed  little  to  the  over-all  program. 
Associated  with  the  conduct  of  the  tests  described  above,  another 
phenomenon  was  made  apparent.  At  some  points  on  the  shock 
wave  where  tests  were  being  conducted  through  large  ranges  of 
rod  angular  variation,  it  was  noted  that  for  large  negative  angles 
of  attack  and  after  boundary  layer  separation  on  the  rod  had  been 
initiated,  the  entire  shock  wave  separated  to  the  forward  end  of 
the  rod.  An  explanation  of  this  phenomenon  is  not  possible  at 
this  time  since  it  does  not  appear  to  be  associated  with  a  center- 
line  disturbance  ahead  of  the  shock  generator,  since  it  occurred 
at  various  rod  end  locations  relative  to  the  centerline  of  the 
shock  generator. 

These  tests  were  immediately  followed  by  the  first  para¬ 
chute  test  program. 

3.  Parachute  Test  Program 

a.  Wind  Tunnel  Test  Installation  and  Deployment  Method 

All  test  models  in  the  Phase  I  program  were  tested  behind 
a  bi-conic  forebody  (Figure  II-5)  which  was  strut  mounted  to  the 
tunnel  side  wall.  A  spring  actuated  deployment  mechanism  was 


ASD-TDR-62-«44 


22 


Figure  II-5  Phase  I  Forebody,  Parachute  Tests 


used  to  deploy  the  models  when  tile  desired  'utmel  flow  conditions 
were  attained.  'I'lie  deployment  mecliainsm  seas  actuated  by  an 
external  meclianic.il  pull-c.iblc-  control.  Provisions  svere  tilso 
made  for  tlie  possible  d-ployiin  iit  of  two  models  during  a  tunnel 
run.  One  model  (witlumt  provision  for  drag  measurements)  ssas 
to  be  stowed  in  a  rip-cover  on  tlie  trailing  edge  of  tiie  strut,  and 
mounted  to  the  forebody  base  by  means  of  a  yoke  and  ball  re¬ 
lease  connection.  An  exti-rn.il  pull  i  o  rd  seas  rigged  so  as  to  run 
from  the  rip-cover  to  outside  the  tunnel.  By  pulling  this  cord, 
tlie  first  parachute  would  be  deployed  ujion  removing  the  rip- 
c  o  ve  r . 

Prosnsion  ssas  m.ide  for  the  first  par.ichute  to  be  jettisoned 
dossTistream  .it  the  conclusion  of  the  test  by  .1  release  mechanism 
triggered  by  an  initial  displacement  of  the  cable;  further  dis¬ 
placement  of  tile  cable  actuated  the  deployment  of  the  second 
parachute  model  stowed  in  the  fore  body.  Due  to  possible  com¬ 
plications  to  the  tunnel  operation  caused  by  a  jettisoned  para¬ 
chute,  only  the  spring  actuated  deployment  system  was  used. 

Drag  measurements  were  provided  by  a  strain  gauge  mounted 
drag  link  installed  inside  the  forebody  whose  circuit  was  connect¬ 
ed  to  the  tunnel  facility  instrumentation  equipment.  Force  meas- 
surements  were  recorded  on  Visicorder  tape.  Figure  II-6  is  a 
photograph  of  the  parachute  lest  installation  used  in  the  Phase  I 
tests. 

b.  Test  Model  Descriptions 

In  the  first  series  of  parachute  tests,  a  study  was  conduc-^ 
ted  involving  five  modified  guide  surface  type  canopies,  five 
conical  type  canopies,  a  hemisflo  canopy,  and  a  flat  circular 
ribbon  canopy  with  the  capability  of  remote  variation  of  reefing 
area  ratio.  This  selection  conforms  with  the  trends  indicated 
in  Section  I-D  and  with  contractual  requirements.  Variations 
in  the  basic  gore  designs  of  the  guide  surface  types  were 
accomplished  in  order  to  provide  different  inlet  to  maximum 
projected  diameter  ratios  as  well  as  inlet  angles,  the  latter  being 
based  upon  analyses  similar  to  that  discussed  in  Section  I-D. 
Specifications  for  the  various  configurations  tested  in  this  pro¬ 
gram  are  shown  in  Table  II-3.  Sketches  and  material  specifi¬ 
cations  of  these  configurations  are  shown  in  Figures  II-7,  11-8, 
and  II-9,  and  Table  II-4, 
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igure  II-6.  Phase  I  Parachute  Test  Installation 


PHASE  I  PARACHUTE  TYPES  AND  SPECIFICATIONS 
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igure  II-7.  Phase  I  Bi-Conic  Parachute  Test  Models 


TABLE  II-4 


PHASE  I  MODIFIED  GUIDE  SURFACE  MODEL  SPECIFIC  AT  IONS 


TYPE 

A&B 

De 

0.  842 

Dj 

5.  89 

Dp 

8.  42 

0 

23.  7° 

1 

8 

L 

3.  51 

a 

0.732 

/3 

oc 

o 

• 

o 

0.  H-12 

4.  21 
8.  42 

22.  2“ 

8 

5.  85 
0.  529 

79.  5° 


0.  842 
2.  95 
8.  42 
21.7° 

8 

7.6 
0.  367 
79.  5° 


0.  842 
4.  21 
8.  42 
12.  2° 

8 

10.  5 
0.  524 
67° 


0.  57 

0.  24 

0.  81 

0.  44 

1.13 

0.71 

1,49 

0.  97 

1. 85 

1.  30 

2.  23 

1.  68 

2.  62 

2.  15 

3.  00 

2.  59 

3.  23 

1.90  3,  46 
2.  38  3.  73 
2.81  3.91 


6.  68 
7.48 
9.92 


ALL  CHUTES 
ARE  MADE 
FROM  FOLLOW- 
INC  MATERIALS; 


GORE:  CLOTH, 
NYLON 
MIL-C-8021 
TYPE  I 


SUSPENSION 
LINE:  CORD, 
NYLON 
MIL-C-5040 
TYPE  I 

RISER:  WEBBING 
NYLON 
MIL-W-4088 
TYPE  I 
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c,  Rir. suits  of  Phase  I  Tests 

In  the  conduction  of  this  testinj*  program,  the  basic  instru¬ 
mentation  and  equipment  employed  in  the  tests  provided  force 
measurements,  direct  motion  photography,  and  Schlicren  movies. 
Analysis  of  the  force  data  recorded  on  Visicorder  tape,  as  well 
as  the  motion  pictures,  resulted  in  both  favorable  and  unfavorable 
performance  characteristics  of  the  test  configurations.  Since  the 
direct  movies  (at  500  frames  per  second)  and  the  Schlieren 
movies  (at  1000  frames  per  second)  made  possible  the  analysis 
of  various  phenomena  at  short  frame-time  intervals,  changes  in 
inlet  opening,  amount  of  canopy  inflation,  steadiness  of  shock 
geometry  and  canopy  stability  characteristics  were  readily  dis- 
cernable.  On  the  basis  of  this  analysis  Table  11-5  was  prepared 
which  summarizes  the  essential  performance  parameters  and 
characteristics  of  the  test  models. 

In  the  accomplishment  of  this  test  program,  indications  of 
favorable  canopy  geometry  were  noted.  In  one  instance  (modified 
guide  surface  type  D)  an  over-ail  indication  of  good  performance 
prevailed,  although  the  drag  was  low  compared  to  types  A  and  B. 
Since  the  suspension  lines  twisted  due  to  model  rotation  within 
the  first  0.  5  second  after  deployment,  a  thorough  study  of  this 
test  canopy  was  not  possible.  During  this  time  interval,  however, 
direct  photography  indicated  excellent  inflation  characteristics 
as  seen  in  Figure  11-10.  Although  the  drag  coefficient  of  0,  095 
based  upon  constructed  cloth  area  at  a  Mach  number  of  Z.  99,  was 
not  high,  the  important  result  of  the  Phase  I  test  program  was 
the  attainment  of  favorable  behavior  of  the  canopy  such  as  infla¬ 
tion  and  stability  characteristics,  and  the  steadiness  of  the  shock 
wave.  The  latter  w'as  not  observable  in  the  case  of  configuration 
D  since  the  line  twist  occurred  prior  to  initiation  of  the  Schlieren 
photography. 

From  the  standpoint  of  drag,  the  modified  guide  surface 
type  A  canopy  appeared  best.  This  test  canopy  provided  a  drag 
coefficient  of  1.  Z4  based  upon  the  design  projected  frontal  area 
or  0.  37  based  upon  the  constructed  cloth  area  at  a  Mach  number 
of  1,  99.  This  model  was  found  to  be  somewhat  unstable  about  its 
center  of  gravity. 

Favorable  characteristics  of  other  canopies  were  observed 
in  close  examination  of  the  high  speed  motion  pictures.  Although 
in  most  cases,  other  than  the  types  A,  B,  and  D,  over-Eill  canopy 
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Notes:  (1)  Suspension  lines  are  twisted 

(2)  Cj-j  based  on  average  maj^imum  inflated  frontal  area 

(3)  Constructed  d  =  8  inches 


Notes;  (1)  Suspension  lines  are  twisted. 

(2)  Cj)  based  on  average  maximum  inflated  frontal  area. 

(3)  Constructed  d^  =  8  inches. 
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Notes:  (1)  Constructed  d^^  =  8  inches 

(2)  Cj-)  based  on  average  majcimum  inflated  frontal  area 


Notes:  (2)  based  on  avera^'e  maximum  inflated  frontal  ar 

(3)  Constructed  d  =8  inches 
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Notes:  (2)  Cj)  baaed  on  average  maximum  inflated  frontal  area 

(3)  Constructed  d  =  8  inches 

(4)  (d.  /2/3dQ)2 


Figure  II-IO.  Phase  I  Modified  Guide  Surface  Model  D.  Film  Speed 


500  fps,  M 


=  2.  99,  q  =  1 56  psf,  x/d 
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behavior  was  not  outstanding.  However,  in  several  configura¬ 
tions  there  were  indications  of  good  inlet  inflation,  roof  per¬ 
formance  and  shock  wave  steadiness.  Outstanding  among  these 
is  the  type  G  conical  canopy.  The  inlet  geometry  for  this  con¬ 
figuration  may  be  seen  (Fig.  11-7)  to  be  similar  to  the  Type  D 
guide  surface  model.  The  inlet  of  this  model  appeared  to  be 
very  well  inflated  during  its  test  when  at  near  zero  angle  of 
attack.  The  Type  F  conical  appeared  best  of  ail  conical  types 
from  the  standpoint  of  stability,  although  the  canopy  was  under¬ 
inflated.  The  conical  type  parachutes  indi'*ated  (except  for  Type 
F)  varying  degrees  of  instability  about  the  point  of  suspension. 
The  reefed  flat  circular  ribbon  canopy  exhibited  very  good  sta¬ 
bility  and  good  roof  inflation,  particularly  when  reefed  to  Z8 
percent  at  a  Mach  number  of  Z,  99. 

In  Table  11-6  are  listed  the  essential  parachute  components 
and  associated  geometry  which  indicated  good  component  or 
over-all  performance  in  the  Phase  I  tests.  The  configurations 
with  which  these  components  or  geometric  properties  are  asso¬ 
ciated  are  also  tabulated. 


TABLE  II-6 


DEMONSTRATION  OF  FAVORABLE  COMPONENT  PERFORMANCE 

PHASE  I  TEST  MODELS 


Component 

Performance 

Inlet  Inflation 

Roof  Inflation 

Drag 

Stability  about  point 
of  suspension 

Stability  about  C.  G. 

’i-'Reefed  circular  flat  ribbon 
^’'■'Hemisflo 


Associated  Test  Model 
Configuration _ 

A,  B,  C,  D,  E,  G 

A,  B,  D,  RR’’' 

A,  B 

A,  B,  C,D,E,F,  RR=f<Hem.** 

C,D,E,F,  RR*,  Hem** 
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It  may  be  noted  also  that  the  number  of  appearances  which 
the  various  configurational  types  made  in  this  tabulation  directly 
indicate  the  number  of  favorable  performance  characteristics 
associated  with  a  particular  model  configuration.  Configura¬ 
tions  A,  B,  and  D  each  appear  four  of  five  times  in  this  tabula¬ 
tion.  Both  A  and  B  fail  to  appear  in  the  category  of  good  stabil¬ 
ity  about  the  model  center  of  gravity.  The  Type  D  configuration 
does  not  appear  in  the  outstanding  drag  classification,  which  is 
expected  since  the  inlet  area  of  this  model  is  much  smaller  than 
that  for  Types  A  and  B,  which  exhibit  very  high  drag. 

The  above  mentioned  results  were  utilized  in  selection  of 
configurations  for  Phase  la  tests  which  are  now  described. 

B.  Phase  la  Parachute  Test  Program 


1.  General  Configurational  Considerations 

This  test  program  intervened  between  the  first  and  second 
Government  facility  testing  programs  and  was  conducted  in  the  Cook 
Technological  Center  wind  tunnel  facility.  The  objective  of  this 
program  was  to  improve  the  Phase  I  configurations.  The  approach 
was  to  evaluate  configurations  which  had  geometric  characteristics 
similar  to  those  models  which  performed  best  in  the  Phase  I  program. 
As  discussed  in  Part  A,  certain  configurations  tested  in  the  Phase  I 
program  indicated  either  favorable  inlet  inflation  characteristics,  or 
over -all  stability,  or  canopy  roof  inflation,  or  ove  r-all  inflation,  or 
over-all  inflation  and  stability.  Since  several  configurations  indica¬ 
ted  favorable  inlet  inflation  characteristics,  these  geometries  were 
retained  in  future  model  designs,  leaving  the  canopy  roof  study  as  the 
primary  emphasis  of  this  program.  The  non-porous  bi-conical  con¬ 
figurations,  having  indicated  under-inflation  characteristics,  particu¬ 
larly  in  the  roof  regions,  were  modified  in  such  a  way  as  to  decrease 
the  exit  to  inlet  area  ratio.  Since  the  type  F  conical  canopy  indicated 
good  stability,  and  also  represented  the  conical  type  with  the  smallest 
total  surface  area,  it  was  considered  for  further  study  in  this  inter¬ 
im  test  program.  Modifications  were  made  in  the  roof  geometry  to 
improve  the  inflation  characteristics. 

The  tests  in  the  Phase  I  program  involving  the  modified  guide 
surface  type  canopies  indicated  excellent  stability  with  reference  to 
the  point  of  suspension,  and  except  for  types  A  and  B,  they  had 
excellent  stability  with  reference  to  their  center  of  gravity.  Inflation 
characteristics  were  fair  to  good  for  some  configurations  (Types  A, 
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B,  C  and  D),  and  as  previously  mentioned,  were  best  for  the  Type 
D.  In  some  instances,  involvinj*  tests  of  these  configurations,  some 
oscillation  of  the  guide  surface  ^>anels  were  noted.  The  tendency  to 
rotate,  resulting  in  suspension  line  wrap-up,  was  also  indicated  in 
those  tests  in  which  a  sv/ivel  was  not  used.  It  was  intuitively  appar¬ 
ent  that  mal-alignment s  of  the  guide  surface  panels  due  to  fabrication 
tolerances  would  tend  to  cause  rotation.  Tests  indicated  this  to  be 
true,  particularly  when  the  guide  surface  panels  were  well  inflated 
(e.g. ,  Type  D).  It  was  not  apparent  from  the  tests  that  the  guide 
surface  construction  contributed  significantly  to  the  favorable  per¬ 
formance  of  the  guide  surface  models.  Accordingly,  part  of  the 
study  in  the  interim  program  was  aimed  at  evaluating  the  contribution 
of  the  guide  surface  panels.  This  was  accomplished  by  constructing 
models  with  geometry  similar  to  a  particular  guide  surface  model. 
Simultaneously,  a  study  of  canopy  roof  geometry  was  conducted. 

The  instability  about  the  center  of  gravity  which  was  character¬ 
istic  of  Types  A  and  B  was  attributed  to  excess  negative  lift  contribu¬ 
ted  by  the  flat  roof.  These  configurations  indicated  excellent  drag  and 
exhibited  good  inflation  and  were  also  constructed  with  the  smallest 
cloth  area  of  the  guide  surface  types  (all  having  the  same  frontal  area). 
Hence,  they  were  considered  particularly  worthy  of  further  study  with 
modifications  to  improve  their  stability  characteristics. 

In  the  tests  of  the  flat-circular  ribbon  model  in  the  Phase  I  pro¬ 
gram,  good  inflation  of  the  canopy  roof  was  indicated  when  the  model 
was  in  a  reefed  state.  Flagging  of  the  forward  ribbons  was  apparent, 
and  expected  due  to  excess  material  in  that  area.  Attention  was  given 
to  the  design  of  a  model  with  a  ribbon  roof,  and  a  shaped  inlet.  The 
inlet  design  was  aimed  at  simulating  the  ribbon  inlet  in  the  best  reefed 
condition  (Z8%  reefing  ratio)  of  the  Phase  I  tests. 

As  discussed  above,  the  primary  emphasis  of  the  study  perform¬ 
ed  in  the  interim  test  program  was  the  application  of  favorable  canopy 
geometry  indicated  by  the  Phase  I  test  models.  The  interim  test  pro¬ 
gram  was  also  to  furnish  data  from  which  an  evaluation  may  be  made 
of  the  contribution  to  parachute  performance  of  the  guide  surface  panel 
construction. 

Z,  Description  of  the  Test  Model  Configurations 

The  models  for  the  interim  test  program  were  designed  in 
accordance  with  the  conclusions  drawn  from  the  performance  analysis 
of  the  Phase  I  test  configurations  and  are  described  in  Table  II-7.  The 
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TABLE  II-7 

PHASE  la  I'EST  MODEL  SPEC  I  PIC  A'lTONS 
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listing  in  the  lablo  c or ros ponds  geni-rally  to  tin-  order  in  which  tlu; 
models  were  tested.  It  is  noted  that  the  configurations  C-1,  C-3, 
C-2,  C-‘l  and  C-S  are  niodili<  ations  of  the  Phase  I  model,  Type  F, 
Model  G.S.  -1  is  a  sealed  duplu  au-  of  the  Phase  I  test  model  Type  A, 
Since  model  Types  A  and  B  are  quite  similar,  only  one  of  these  types 
(A)  was  employed,  along  with  its  modifications  in  the  Phase  la  study. 
Test  models  G.S.  -Z,  G.S.  --1  and  G.S.  -5  are  modifications  (in  roof 
configurat  ioti)  of  model  G.  S.  - 1 .  'fhe  test  model  C.R.-l  is  a  simula¬ 
tion  of  the  reefed  ribbon  configuration  which  appeared  to  perform 
quite  well  in  the  Phase  I  tests.  Test  model  C-6a  was  designed  to 
duplicate  the  general  g«-ometry  of  model  G.  S.  -5  with  the  purpose  of 
evaluating  the  e ffet  t i vene s s  of  the  guide  surface  panel  construction 
employed  in  tlie  model  G.S.  -  S.  The  test  model  C-9  w’as  designed  to 
be  tile  conical  equivalent  of  the  Phase  I  i»  st  morlel  Type  D  for  a 
similar  evaluation.  Models  C-7a  and  C-7b  were  designed  for  further 
study  of  roof  geometry  and  i-xit  tt>  inlet  area  ratio.  These  models  are 
further  modifications  of  test  model  C -6a.  Photographs  of  test  models 
C-7a  and  C-6a  are  shown  in  Figure  II-ll.  Details  of  all  Phase  la 
models  are  given  in  Figures  II-H,  11-15,  II-16,  11-17,  11-18,  and 
Tables  11-8  and  II-9. 

3.  Wind  Tunnel  Test  Installation  and  Deployment  Methods 

In  the  interim  test  program  in  the  Cook  Technological  Center's 
wind  tunnel  facility,  all  parachute  models  were  tested  in  the  wake  of 
a  cone-cylinder  forebody  which  was  mounted  in  the  test  section  as 
shown  in  Figure  11-12.  Test  models  employed  risers  which  located 
the  model  skirt  approximately  ten  forebody  base  diameters  down¬ 
stream.  High  speed  (1  000  fps)  Schlieren  photography  data  was  ob¬ 
tained  during  each  test.  Force  data  were  not  obtained  during  this 
series  of  tests  since  the  primary  objective  of  the  study  was  the  eval¬ 
uation  of  the  behavior  of  the  model  from  the  standpoints  of  inflation, 
stability  about  the  point  of  suspension,  and  characteristics  of  the 
shock  geometry  ahead  of  the  model.  The  ratio  of  model  diameter  to 
forebody  diameter  (2,  18)  employed  in  these  tests  approximate  that 
used  in  the  Phase  1  tests  (2.  13). 

Parachute  models  were  stowed  in  a  bag  secured  by  a  break 
thread  to  the  forward  end  of  the  riser  (near  the  forebody  base)  during 
tunnel  run-up.  When  flow  conditions  in  the  tunnel  were  established, 
the  bag  was  pulled  off  by  means  of  a  line  running  from  the  bag  to  a 
downstream  mechanism  capable  of  axial  translation  by  remote  manual 
control. 
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Figure  II-ll.  Phase  la  Test  Models  C-6a  and  C 
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II  12  Phase  I  a  Porebody,  Parachute  Tests 


igure  11-14.  Phase  la  Bi-Conic  Parachute  Models 
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This  system  is  illustr.ited  in  Figure  11-13,  When  sufficient  tension  in 
the  pull  line  was  applied,  failure  of  the  parachute  deployment  bag  break¬ 
line  resulted  sn  deployment  of  the  model.  The  mechanism  used  to  deploy 
the  model  is  part  of  the  installation  used  in  the  Cook  Technological 
Center's  Basic  Test  Program,  discussed  in  Section  1. 

4.  Results  of  the  Test  Program  and  Data  Analysis 

In  the  conduct  of  this  test  series,  initial  model  designs  were  per¬ 
formed  in  accordance  '.^.ith  general  modifications  as  indicated  by  the 
analysis  of  the  Phase  1  test  program  and  discussed  previously.  However, 
structural  and  fabrication  problems  developed  as  the  test  series  pro¬ 
gressed,  requi ring  construction  modifications  as  indicated  by  test 
model  performance  analysis,  and  as  required  for  test  model  survival 
throughout  the  testing  period.  The  high  dynamic  pressure  (700  psf) 
environment  of  these  tests  imposed  quite  severe  requirements  on  the 
materials  employed,  as  well  as  the  fabrication  techniques.  However, 
minimum,  weight,  strength  and  stiffness  materials  were  employed  where 
possible  in  thu  fabric.ition  of  the  models  in  order  to  maintain  as  limp 
a  construction  as  possible.  In  all  tests,  the  model  canopies  were  de¬ 
ployed  from  a  tightly  packed  initial  condition,  requiring  that  they  self- 
inflate  after  Ueploymi  nt. 

In  Tabic  !1-1U  the  tests  conducted  in  this  series  are  listed  chron¬ 
ologically.  Siabil'.ty  about  the  point  of  suspension  (Aa°)  and  inflation 
characteristic.s  are  lab.ilated.  It  may  be  noted  that  successive  tests 
do  not  necessarily  fully  explore  a  given  configuration  along  with  modi¬ 
fications  indicated  by  the  previous  test.  This  occurs  because  study  of 
the  data  derived  from  a  given  test  was  not  conducted  before  testing  of 
a  modification  of  this  canopy.  During  this  study  period,  different  con¬ 
figurational  types  were  tested. 

Tests  of  canopies  which  were  constructed  of  non-porous  materials 
were  found,  as  in  the  Phase  I  tests,  to  suffer  structural  damage  either 
early  in  the  test  period  or  before  wind  tunnel  shutdown.  Model  canopies 
in  the  Phase  I  tests  as  well  as  models  C-1,  C-Z,  C-3,  and  C-4  in  the 
Phase  la  tests  were  constructed  of  Dantex  ,  a  neoprene  coated  nylon 
fabric,  A  heavier  material,  Vulcan,  was  also  considered  for  the  con¬ 
struction  of  the  non-porous  canopies,  and  was  employed  for  one  of  the 
later  (C-5)  tests.  However,  as  indicated  previously,  it  was  desired 
to  construct  canopies  in  as  limp  a  form  as  possible,  and  resort  to  this 
heavier  material  was  not  made  until  conclusive  evidence  of  its  need 
was  shown  by  tests.  It  may  be  appropriate  to  point  out  at  this  time 
that  static  loadings  are  not  critical  since  even  the  use  of  lightweight 
materials  in  model  canopy  construction  result  in  high  static  factors  of 
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TAP.LE  II- 10 


PHASE  la  PARACHUTE  TEST  SUMMARY 


onfigu  ration 

oo 

Inflation 

C  -  la 

2.  10 

2.  5 

Good 

G.S.  -1 

2.  OH 

15.  0 

Good 

C.  R.  -1 

2.  1  1 

0.  0 

Poor 

G.S.  -2 

2.  12 

8.  0 

Good 

C-4 

Subs 

- 

Good  subsonic 

G.S.  -4 

2.  09 

10.  0 

Good 

G.S.  -5 

2.  10 

9.  0 

Good 

C-5 

2.  01 

7.  0 

Good 

C  -6a 

2.  1  1 

5.  0 

Good 

C-9 

2.  06 

1.  5 

F  air 

C-7b 

2.  06 

3.  0 

Good  to  Poor 

C-7a 

2.  09 

5.  0 

Good  to  Poor 

safety.  Instead,  dynamic  loadings  caused  by  high  frequency  oscilla¬ 
tions  appear  to  be  critical,  and  since  these  are  difficult  to  predict, 
material  requirements  must  be  based  upon  test  results. 

The  first  test  in  the  Phase  la  series  involved  the  model  canopy 
C-1  which  was  designed  to  consist  of  a  roof  modification  of  the  Phasel 
Type  F  conical.  The  roof  angle  was  increased  (to  keep  the  total  con¬ 
structed  area  low)  and  the  vent  area  decreased  in  order  to  decrease 
the  exit  to  inlet  area  ratio.  The  reduction  in  exit  to  inlet  area  ratio 
was  aimed  at  improving  the  inflation  characteristics  of  the  Type  F 
canopy  which  indicated  under -inflation  in  the  Phase  I  tests.  Although 
the  C-1  model  deployed  prematurely,  it  indicsited  good  inflation  in  the 
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subsonic  flow  conditions  prfvailin{>  durinj>  tunnel  shutdown.  Examina¬ 
tion  of  the  canopy  after  the  test  revealed  initiation  of  structural  failure, 
and  it  was  decided  that  the  canopy  with  tliis  low  (0.  16)  exit  to  inlet  area 
ratio  would  not  survive  supersonic  test. 

The  model  C-  ia  test  followed.  With  an  exit  to  inlet  area  ratio 
of  0.  50,  very  good  inflation  was  indicated  (Figure  11-19),  although  a 
Very  low  frecpieiu  y  breathing  was  evident.  The  model  was  highly  stable, 
with  oscillation  limits  about  the  j)oint  of  suspension  of  ±2.  5°.  The  low 
frequency  inflation  and  collapse  of  tins  model  is  considered  to  indicate 
marginal  inflation  characteristi<  s  for  the  Mach  number  (M  =  2.  1 )  of 
this  test.  The  canopy  typo  C-2  ti-st  indicated  fair  stability  and  good 
inflation  although  examination  of  the  Schlicren  film  indicated  canopy 
damage  shortly  after  deployment.  This  canopy  had  approximately  the 
same  exit  to  inlet  area  ratio  as  the  type  C- !>,  but  with  high  inlet  to 
maximum  diameter  ratio  (0.  9  compared  to  0.  7). 

Test  model  G.S.  -1  was  designed  to  duplicate  the  Phase  I  test 
model  modified  guide  surface  Type  A,  which  indicated  excellent  infla¬ 
tion  in  the  Phase  I  tests.  Since  modifications  of  the  roof  geometry 
were  contemplated  in  order  to  improve  stability  characteristics  about 
the  canopy  center  of  gravity,  this  configuration  was  scaled  to  the 
Phase  la  model  size  and  tested  for  comparative  purposes.  As  in  the 
Phase  I  tests,  the  canopy  inflation  proved  to  be  excellent  with  instability 
about  the  point  of  suspension  of  ±15°.  Instability  about  the  canopy  cen¬ 
ter  of  gravity  was  not  apparent  in  this  test,  due  possibly  to  differences 
m  canopy  suspension  system  mass  distributions  as  compared  to  the 
Phase  I  test. 

Test  model  C.R. -1  was  constructed  with  a  non-porous  inlet,  and 
a  hemispherical  ribbon  roof.  The  gore  patterns  for  the  roof  were  con¬ 
structed  from  a  16-gore,  22.  5  percent  geometric  porosity  hemisflo 
design.  This  configurational  study  was  suggested  by  the  Phase  I  test 
performance  of  the  reefed  circular  flat  ribbon  canopies,  where  good 
inflation  and  stability  were  shown  at  M  =  3.  The  shaped  conical  inlet 
was  employed  to  eliminate  the  excess  skirt  region  material  normally 
associated  with  a  reefed  parachute.  In  the  test  of  this  model,  very 
poor  inflation  and  severe  breathing  were  evidenced.  The  poor  inflation 
may  be  attributed  to  the  high  (0.  80)  exit  to  inlet  area  ratio. 

Model  G.S.  -2  was  designed  to  provide  a  reduction  in  the  flat  plate 
lift  associated  with  model  G.S.  -1,  by  employing  an  annular  vent.  This 
canopy  proved  to  be  very  well  inflated,  and  exhibited  a  reduction  in 
instability  about  the  point  of  suspension  (  8°)  as  compared  to  model 
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Fieure  11-19.  Schlieren  Photograph  of  Phase  la  Model  C-3a 
Film  Speed:  l.OOOfps,  M  =  2.10,  q  =  705ps£, 
x/d  =  10 


I 

Figure  11-20.  Schlieren  Photograph  of  Phase  la  Model  C-6a.  Film  Speed: 

1, 000  fps,  M  =  2. 1  3,  q  =  695  psf,  x/d  =  10 
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G.S.  -1  (15°).  Analysis  of  the  Schlicron  films  for  this  test  indicated 
a  considerable  amount  of  relative  motion  of  the  roof  center  disc  with 
respect  to  the  main  canopy.  Relative  motion  of  this  type  couJd  result 
in  an  unsymiT.etrical  vent  area  distribution,  hence,  a  further  reduction 
in  canopy  oscillation  about  the  point  of  suspension  should  residt  by 
providing  a  fixid  vent  area.  It  was  decided  to  retest  this  configuration, 
giving  some  attention  to  improving  the  roof  disc  support. 

The  test  model  C-9  (roof  modification  of  Type  F)  was  studied 
next  in  the  test  series.  This  model  has  a  similar  vent  construction  to 
the  model  G.S.  -L.  -\lthough  premature  deployment  resulted  in  this 
test,  the  instability  of  the  roof  disc  was  marked,  as  in  the  case  of 
test  model  G.S.  -2.  The  use  of  a  porous  mesh  material  or  a  geodetic 
type  of  line  suspension  from  the  canopy  to  the  disc  to  pro\dde  a  uni- 
forn  exit  geometry  was  considered  in  succeeding  canopy  designs 
employing  an  annular  vent. 

In  order  to  test  the  capability  of  a  non  "mil-spec"  mesh  material 
and  also,  the  merits  of  a  central  vent,  the  test  model  G.S.  -4  was 
designed.  In  this  model,  a  40  percent  porosity  nylon  mesh  was  used. 
The  mesh  material  filled  nearly  the  entire  roof,  providing  a  central 
vent.  In  the  test  of  this  model,  very  good  inflation  was  realized,  and 
oscillation  angles  about  the  point  of  suspension  (10°)  indicated  better 
stability  than  the  G.  S.  -1  (1  5°)  but  poorer  than  the  model  G.  S.  -  2  (8°). 
The  exit  to  inh-t  .irea  ratio  of  this  configuration  is  approximately  the 
same  as  that  for  the  model  G.S.  -2, 

Canopv  model  G.S.  -5  was  designed  to  have  the  same  exit  to  inlet 
area  ratio  -s  model  G.S.  -4,  but  with  an  annular  vent  filled  with  40 
percent  porosity  mesh  to  secure  the  roof  center  disc  and  maintain  a 
fixed  vent  geometry.  Analysis  of  the  Schlieren  films  indicated  that  a 
reduced  oscillation  (9°)  about  the  point  of  suspension  resulted  with  this 
configuration  as  compared  to  canopy  models  G.S.  -1  and  G.S.  -4.  Ex¬ 
cellent  inflation  of  the  canopy  was  realized,  and  the  use  of  the  mesh 
material  appeared  to  provide  a  more  uniform  vent  geometry. 

Configuration  C-5  was  designed  to  represent  a  modification  of 
the  type  F  conical  model.  A  flat  roof  with  an  annular  vent  providing 
an  exit  to  inlet  area  ratio  of  0.  50  was  used.  This  model  was  con¬ 
structed  of  non-porous  Vulcan  material.  A  different  means  of  secur¬ 
ing  the  roof  center  disc  was  attempted  in  this  model  by  employing  a 
bicycle  spoke  arrangement  of  lines  from  the  canopy  to  the  disc. 
Schlieren  film  analysis  indicated  small  (7°)  canopy  oscillations  about 
the  point  of  suspension.  Rapid  breathing  and  considerable  relative 
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motion  of  tlif  roof  dirfo  with  rcspfct  to  the  m.iin  canopy  were  obaerved. 
The  disc  motion  reanltfd  in  attachment  line  failure  and  ultimate  de¬ 
struction  of  the  mode!  [)rior  to  tunnel  shutdown. 

Model  canopy  C-ba  was  designed  in  order  to  evaluate  the 
relative  performances  of  a  conical  type  canopy  with  the  guide  surface 
type  construction.  Tlur  model  was  designed  to  duplicate  the  geometry 
of  the  model  G.S.  -S  in  inlet  to  maximum  diameter  ratio,  inlet  angle, 
roof  geometry,  exit  to  inlet  area  ratio  and  construction  material.  The 
model  was  constructed  with  a  lower  exit  to  inlet  area  ratio  than  that 
of  G.  S.  -5,  however.  A  much  lower  oscillation  angle  (5°)  residted  in 
the  test,  with  excellent  inflation  (Figure  11-20).  This  appeared  to 
indicate  that  there  is  no  specific  contribution  of  the  guide  surface  con¬ 
struction,  except  p<-rhaps  the  unfavorable  rotation  tendency  of  this 
type  of  construction  c  aused  by  asymmetries  introduced  in  fabrication. 

The;  model  canopy  C-9  was  next  tested  in  the  Phase  la  series. 
This  canopy  was  designed  to  b<‘  the  conical  equivalent  of  the  design 
geometry  of  the  Phase  I  modified  guide  surface  model  D.  The  Type 
D  model  was  described  in  Section  1,  and  as  indicated  in  that  section, 
proved  to  be  the  best  performing  configuration  in  the  Phase  I  tests 
from  the  standpoint  of  inflation  and  stability  characteristics.  In  the 
Phase  la  test  of  the  model  C-9,  the  canopy  was  not  well  inflated, 
although  highly  stable  (1.  5°  oscillation). 

It  was  decided  to  design  two  model  canopies  with  moderately 
high  exit  to  inlet  area  ratios  (0.  5  and  0.  6)  to  determine,  if  possible, 
the  limiting  area  ratio  for  good  inflation  at  the  M  =  2  condition.  Test 
models  C-7a  and  C-7b  were  designed  with  this  in  view.  Model  C-7b 
was  designed  with  the  same  inlet  geometry  as  models  G.S.  -1,  G.S.  -2, 
G.S.  -4,  G.S,  -5  and  C-6a,  The  roof  is  flat  and  consists  of  a  35  per¬ 
cent  porosity  mesh  (Perlon)  monofilament.  The  suspensb  n  lines  are 
not  carried  over  the  roof,  as  in  the  case  of  C-7a.  Examination  of  the 
Schlieren  films  of  the  C-7b  test  reveals  a  strong  shock  in  front  of  the 
canopy  when  fully  inflated  with  indications  of  severe  shock  wave  sus¬ 
pension  line  boundary  layer  interaction.  This  type  of  shock  wave 
structure  appears  to  indicate  unstable  flow  conditions  ahead  of  the 
canopy.  This  was  also  observed  in  the  test  of  Model  C-3  previously 
described.  A  very  low  frequency  canopy  inflation  and  collapse 
(breathing)  appeared  to  exist  in  the  C-7b  model  test. 

Schlieren  films  of  the  Model  C-7a  canopy  test  revealed  a  delayed 
opening  of  the  canopy  followed  by  a  period  of  excellent  inflation, 
stability  (5°)  and  shock  wave  structure  (Figure  11-21).  Upon  an  angular 
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Figure  11-21.  Schlieren  Photograph  of  Phase  la  Model  C-7a. 

Film  Speed:  1, 000  fps,  M=  2.  09.  q  =  704  psf,  x/d  =  10 
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excursion  of  approximately  5^^  about  the  point  of  suspension,  which 
brings  the  canopy  in  close  proximity  to  tlie  shock  wave,  the  canopy 
collapsed  and  remained  so  for  the  remainder  of  the  test  film.  It 
appears  tha.  the  inflation  characteristics  of  this  configuration  are 
marginal  (for  M  =  2),  causing  it  to  be  extremely  sensitive  to  small 
disturbances.  In  contr.ist  with  this,  examples  of  strong  inflation 
cliaracte  r  ist  1C  s  are  represented  by  the  model  C-6a  which  remained 
svell  -  inflated  and  w.is  ms  tisitive  to  disturbances  even  up  to  10*^  of 
angular  oscillation,  aiui  model  G.S.  -1  which  remained  well-inflated 
throughout  angular  excursions  of  1 

A  signific  ant  itu  lu  ase  in  ili  -  effective  exit  to  inlet  area  ratio  of 
several  configurations  was  revealed  by  calculating  this  ratio  on  the 
basis  of  the  inflated  geometry  of  the  canopy  as  revealed  by  the 
Seblieren  photographs.  For  example,  the  model  C-7b  was  constructed 
with  an  c’xit  (open)  to  inlet  area  ratio  of  0.  46.  During  test,  calculations 
based  on  surface  area  increase  reveal  the  effective  exit  to  inlet  area 
ratio  to  be  0.  51.  Similarly,  the  model  C-7a  which  was  constructed 
with  an  exit  -  inlet  area  ratio  of  0.  42  had  an  effective  exit  -  inlet  area 
ratio  of  0.  52.  The  high  dynamic  pressure  encountered  during  these 
tests  (the  order  of  700  psf)  w'as  responsible  for  this  material  stretch. 

The  effect  of  exit  -  inlet  area  ratio  on  stability  is  shown  in 
Figure  11-22  where  the  stability  appears  generally  to  vary  directly 
with  this  ratio  for  a  given  Mach  number,  as  indicated  by  early  tests. 
Also  shown  on  this  figure  are  several  effective  exit  -  inlet  area  ratios 
which  are  based  upon  the  inflated  geometry.  An  analytical  treatment 
of  supersonic  parachute  stability  is  beyond  the  scope  of  the  present 
study. 


In  addition  to  the  exploratory  tests  conducted  in  this  program 
which  led  to  the  final  configurational  designs  of  high  performance 
parachutes,  a  parachute  which  represented  a  parallel  program  of 
development  at  the  University  of  Minnesota  under  the  direction  of 
Prof.  Helmut  G.  Heinrich,  was  tested.  This  configuration  (Figure 
11-23)  consists  of  a  solid  cone  to  which  a  large  vent  flexible  canopy  is 
attached.  Successful  tests  of  this  configuration  have  been  conducted 
by  the  University  of  Minnesota  at  Mach  2  and  Mach  3.  Such  tests  have 
been  normally  conducted  without  a  significant  forebody,  and  have  re¬ 
sulted  in  a  stable  shock  configuration  with  an  attached  shock  at  the 
cone  vertex.  The  test  in  the  Cook  Technological  Center  wind  tunnel 
of  a  4  inch  maximum  diameter  model  in  the  wake  (x/d  =  10)  of  a  2.  75 
inch  cone-cylinder  forebody  (Figure  11-12),  provided  a  parachute 
diameter  to  forebody  diameter  ratio  of  1.  45.  In  this  test,  the  canopy 
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Figure  11-22.  Stability  of  Phase  la  Test  Models  vs.  Exit  to 
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igure  11-23.  Phase  la  Spiked  Parachute  Model 


did  not  inflate  properly,  although  it  appeared  that  a  proper  selection 
of  cone  angle  and  cone  position  woiJd  provide  the  necessary  geometry 
for  a  stable  shock  pattern  if  the  wake  flow  properties  for  the  x/d 
position  at  which  the  cone  and  canopy  are  positioned  are  considered. 
The  primary  advantage  of  this  configuration  consists  of  the  application 
of  a  parachute -like  decelerator  of  large  diameter  to  forebody  diameter 
ratio  at  supersonic  speeds.  As  will  be  discussed  later,  the  favorable 
effects  of  a  forebody  wake  are  essential  to  the  satisfactory  perform¬ 
ance  of  a  trailing  parachute  decelerator  at  supersonic  speeds.  The  use 
of  the  cone  is  intended  to  provide  the  necessary  favorable  flow  condi¬ 
tions  ahead  of  the  canopy  for  satisfactory  operation.  It  may  be  noted 
(Figure  11-23)  that  a  distinct  similarity  exists  between  this  canopy 
inlet  shape  and  that  of  configurations  C-6a  and  C-7a. 

For  parachute  configurations  which  are  constructed  without  a 
conical  appendage  such  as  is  used  in  the  "spiked  parachute  or  cone- 
cup"  configuration  described  above,  satisfactory  wake  conditions  may 
be  provided,  for  example,  by  means  of  a  small  parachute  in  the  wake 
of  which  a  larger  parachute  is  deployed. 

Analysis  of  the  Phase  la  tests  indicated  that  the  general  use  of  a 
low  porosity  material  in  the  canopy  inlet  region  provides  the  necessary 
high  inflation  characteristics  if  the  maximum  inlet  angle  which  is 
employed  is  limited  to  that  which  corresponds  to  the  local  flow  direc¬ 
tion  ahead  of  the  canopy.  The  length  of  the  stabilizing  inlet  conical 
surface  neeo  not  be  large  to  provide  stability  for  the  low  total  canopy 
porosities  tested.  The  conical  extension  used,  however,  is  greater 
(even  for  models  C-6  and  C-7)  than  that  previously  employed  on  con¬ 
ventional  extended  skirt  parachutes.  The  use  of  the  flat  roof  of  high 
relative  porosity  minimizes  the  area  of  the  roof  region  and  results  in 
a  low  ratio  of  total  cloth  area  to  projected  frontal  area.  The  need  for 
shaping  a  canopy  in  the  direction  of  the  expected  inflated  configuration 
appears  unnecessary,  since  configurations  such  as  C-6  and  C-7  appear 
to  adapt  readily  to  a  near-hemispherical  shape  when  inflated. 

The  ratio  of  inlet  diameter  to  maximum  diameter  may  be  large, 
resulting  in  correspondingly  high  drag  coefficients. 

The  contribution  of  porosity  distribution  in  the  roof  area  was 
not  clearly  apparent  as  to  its  effect  on  over-all  performance  at  the 
conclusion  of  the  Phase  la  test  program. 
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c. 


Phase  II  Paraclmte  Test  Pro(»ran) 


I.  Phase  II  Test  Facility 

This  test  phase  was  conducted  in  the  Langley  Research  Center, 
large  supersonic  lunnt;l  facility.  The  test  section  employed  was  the 
4  foot  X  4  foot  Section  L  which  is  capable  of  a  Mach  number  range  of 
1.  3  to  4.65.  A  description  of  this  test  facility  is  given  in  Appendix  UI. 

L,  Selection  of  the  Phase  II  Test  Configurations 

The  results  of  the  performance  analysis  of  the  various  configura¬ 
tions  in  the  Phase  1  and  Phase  la  test  programs  provided  the  basic 
canopy  shapes  which  were  selected  for  the  Phase  II  test  program.  In 
order  to  further  evaluate  the  effects  of  the  forebody  wake  on  the  para¬ 
chute  performance,  model  parachutes  with  inflated  diameters  of  six 
inches  in  some  models  were  selected  to  provide  a  comparison  with  the 
eight  inch  models.  An  inflated  diameter  of  eight  inches  provides  a 
parachute  to  fore  body  diameter  ratio  of  approximately  3,  3  3  and  a  para¬ 
chute  inflated  diameter  of  six  inches  provides  a  2.  5  ratio  of  parachute 
diameter  to  forebody  diameter.  Although  the  model  C-9  did  not  per¬ 
form  to  expectations  in  the  Phase  la  tests,  it  was  decided  to  test  this 
configuration  in  the  Phase  11  program  because  of  its  superior  per¬ 
formance  in  the  Phase  I  tests.  Also,  in  these  latter  tests,  Schlieren 
movies  were  not  obtained  since  the  suspension  lines  of  model  D 
twisted  due  to  rotation  of  the  model  before  Schlieren  coverage  was 
initiated. 

It  was  also  decided  that  a  ribbon  roof  configuration  be  tested  in 
the  Phase  II  program.  The  model  C.R.  -1  had  not  inflated,  as  dis¬ 
cussed  previously,  due  to  the  high  exit  to  inlet  area  ratio.  Further¬ 
more,  the  use  of  the  hemisflo  roof  design  contributed  considerable 
bulk  to  this  portion  of  the  canopy.  In  increasing  the  inlet  area,  by 
using  a  0.  9  inlet  to  maximum  diameter  ratio,  and  constructing  a  flat 
circular  ribbon  roof  of  porosity  16  percent,  where  the  constructed 
diameter  is  the  same  as  the  inflated  canopy  diameter,  an  exit  to  inlet 
area  ratio  of  0.  32  resulted.  Except  for  the  roof  porosity  distribution, 
the  geometric  characteristics  of  this  configuration  are  similar  to  the 
model  C-6a.  This  type  of  roof  construction  would  be  readQy  applica¬ 
ble  to  full  scale  applications  if  found  satisfactory. 

The  model  C-6a,  which  had  superior  performance  in  the  Phase 
la  tests,  along  with  a  modified  (higher  porosity  vent)  model  (0-6  b) 
were  selected  for  the  Phase  II  test  program  in  order  to  evaluate  the 
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performance  of  this  model  up  to  Mach  *1.  Model  C-7a  was  selected, 
due  to  its  marginal  inflation  characteristics  at  Mach  2,  in  order  to 
evaluate  the  effects  on  its  performance  with  increasing  Mach  number. 
This  model  was  adaptable  to  a  reduction  in  exit  to  inlet  area  ratio, 
due  to  its  mesh  roof  construction,  if  the  high  Mach  number  tests 
indicate  a  requirement  :or  this  trend,  althougii  it  was  realized  that 
testing  under  lower  dynamic  pressure  conditions  would  result  in  a 
lower  effective  exit  to  inlet  area  ratio.  By  the  addition  of  a  non-porous 
central  roof  disc  the  maximum  diameter  of  which  is  limited  only  by 
the  diameter  of  the  model,  a  wide  range  of  exit  -  inlet  area  ratios 
could  be  provided.  The  model  C-3  was  selected  because  of  its  good 
performance  in  the  Phase  la  tests,  and  also  because  of  the  good  per¬ 
formance  of  the  conical  Type  F  of  which  it  is  a  modification,  in  the 
Phase  I  tests.  The  model  C-3a  was  further  altered  (C-3b)  for  the 
Phase  II  tests  by  the  addition  of  a  45  percent  porosity  mesh  which 
fills  the  vent,  and  the  vent  diameter  was  increased  so  as  to  provide 
the  original  exit  to  inlet  area  ratio.  Figure  11-24  is  a  photograph 
of  the  Phase  II  test  models. 

The  design  parameters  of  the  Phase  II  test  models  arc  shown  in 
Table  II- 11  and  sketches  and  specifications  of  these  configurations 
are  shown  in  Figure  11-25  and  Table  11-12.  In  addition  to  the  models 
resulting  from  the  Phase  la  tests,  a  conventional  hemisflo  type  para¬ 
chute  was  also  tested.  This  is  the  same  model  used  in  the  Phase  I 
tests,  and  was  retested  to  substantiate  the  results  of  that  program. 

The  model  C -6a  was  tested  employing  both  the  normal  2  dp 
C-6a  (1)  suspension  lines  as  well  as  1  dp  lines  C-6a  (VIII).  This  con¬ 
figuration  was  also  tested  with  both  8  inch  and  6  inch  maximum  diame¬ 
ters.  Test  model  C.  R.  -2  was  tested  employing  both  8  inch  (C.R.  -2b) 
and  6  inch  (C.  R.  -2a)  models.  Configurations  C-3b,  C-6b  and  C-9 
were  tested  employing  only  8  inch  maximum  diameter  models  and  con¬ 
figuration  C-7a  was  tested  by  employing  a  6  inch  maximum  diameter 
model. 

The  canopy  shapes,  along  with  indicated  modifications,  which 
were  considered  for  final  testing  in  this  program  do  not  necessarily 
prescribe  specific  limitations  to  the  geometric  parameters  for  satis¬ 
factory  performance  over  the  range  of  Mach  numbers  considered.  The 
optimum  geometry  for  a  decelerator  application  must  necessarily 
represent  a  compromise  in  drag  effectiveness,  acceptable  stability 
criterir.,  and  system  weight.  In  anticipation  of  attaining  maximum 
drag  effectiveness  with  low  canopy  weight,  the  general  trend  in  the 
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PHASE  II  TEST  MODEL  DESIGN  PARAMETERS 
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igure  11-24.  Phase  II  Test  Models 
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selection  of  lest  inociil  >h.ijjea  for  tlu-  Phase  II  lest  program  was  to 
explore  those  canopies  which  have  higri  inlet  to  maximum  diameter 
ratios,  and  lou  ratios  of  constructed  to  frontal  areas. 

i.  Wind  Tunnel  Installation  atid  Deployment  Method 

Test  models  wer.-  deployed  in  the  wake  of  a  c one  - cylinde r  fore- 
body  (F'igure  II-Z6).  In  order  to  jirovide  tor  the  maximum  use  of  wind 
tunnel  time  during  a  given  tunnel  run-up,  use  was  made  of  a  remotely 
controlled  motor  driven  winch  which,  with  the  double  strut  mounted 
forebody,  was  ecpiijimi  nt  of  the  NASA  Langley  Large  Supe r sonic  Tunnel. 
This  equipment  made  possible  the  rajiid  location  of  the  test  model  at  a 
selected  distance  aft  of  the  forebody  base. 

Before  deploymentjmodcls  were  stowed  in  a  deployment  bag 
which  was  secured  near  the  forebody  base  with  a  break  thread.  The 
bag  was  connected  by  a  line  to  a  downstream  tunnel  strut  which  is 
capable  of  axial  movement,  (Figure  II-Z7).  To  deploy  the  model,  the 
externally  controlled  strut  was  moved  downstream  until  the  break  line 
failed,  allowing  the  parachute  suspension  lines  to  pay  out  and  the  bag 
to  go  downstream,  deploying  the  parachute.  This  same  general  meth¬ 
od  of  deployment  had  proved  successful  in  the  Phase  la  tests  in  the 
Cook  Technological  Center  wind  tunnel  described  in  Part  B. 

4.  The  Phase  II  Parachute  Test  Program 

The  completed  Phase  II  test  schedule  is  shown  in  Table  11-13. 

At  each  test  Mach  number,  data  was  taken  at  either  two  or  three  cano¬ 
py  positions  aft  of  the  forebody  base.  Both  force  measurements  and 
high  speed  (1000  fps)  Schlieren  movies  were  performed  at  each  condi¬ 
tion.  The  listing  in  this  table  is  not  in  the  chronological  order  oftests 
as  performed.  In  general,  models  were  deployed  at  the  high  Mach 
number  which  was  successively  reduced  with  the  model  deployed.  'In 
the  high  Mach  number  range  (M  =  4  to  4.65)  test  models  were  deployed 
at  Mach  4.  Although  the  large  supersonic  tunnel  of  the  Langley  Re¬ 
search  Center  is  capable  of  a  Mach  number  range  of  2.  3  to  4.65, 
three  modes  or  compressor  arrangements  are  necessary  to  accom¬ 
plish  this  range.  After  deploying  and  testing  a  model  in  a  given  mode 
(e.  g. ,  M  =  2.  3  to  2.  75),  rather  than  accepting  the  loss  of  time  asso¬ 
ciated  with  the  tunnel  shutdown  and  redeployment  of  the  model  in  the 
next  higher  mode  (M  =  3.  0  to  3.  75),  the  model  was  left  in  the  tunnel 
requiring  it  to  survive  the  shock  excursion  through  the  test  section 
during  the  mode  change.  This  was  accomplished  successfully  without 
model  damage  in  all  but  two  cases  in  the  entire  test  program,  and 
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TABLE  11-13 


PHASE  II  -  PARACHUTE  TEST  SUMMARY 


Max.  Dia. 

Mach 

9oo 

Configuration 

(In.) 

No. 

X  /d 

(psf) 

C-6a  (I) 

8 

2.  3 

8 

197 

10 

2.  75 

9 

136 

10 

13 

3.  0 

9 

123 

10 

3.2 

9 

103 

10 

3.  5 

9 

81 

10 

3.  75 

8 

65 

10 

4.  0 

8 

159 

10 

C-6b(U) 

8 

2.  75 

9 

152 

10 

3.  00 

9 

123 

10 

3.  20 

9 

105 

10 

3.  5 

9 

81 

10 

3.75 

8 

66 

10 
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TABLE  11-13  (CONT'D) 


PHASE 


Configuration 


C.R.  -2a  (VII) 


C-6a  (V) 


-  PAJ^ACHUTE  TEST  SUMMARY 


Max.Dia.  Mach 


(In.)  No.  x/d 


6  2.  3  7 

10 

12 

2.7  5  7 

10 
12 

3.  0  7 

10 
13 

3.  2  7 

10 
12 

3.  5  7 

10 
12 

3.  75  7 

10 
13 

6  2.  5  7 


10 

1.3 

7 

10 

13 

2. 75  7 

7 

10 

13 


loo 

(£sf) 

214 

152 

123 

105 

81 

66 

220 

117 

96 

158 
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TABLE  11-13  (CONT'D) 


PHASE 

Configuration 
C-6a(V)  (cont’d.  ) 


C-7a  (VI) 


-  PARACHUTE  TEST  SUMMARY 


Max.  Dia.  Mach 


(In.  )  No.  x/d 


6  3.  0  7 

10 

13 

3.  2  7 

10 
13 

3.  5  7 

1  0 
12 

3. 75  7 

10 

13 

10 

4.  0  7 

10 

4.  3  7 

9 

12 

4.65  7 

9 

12 

6  2.  3  7 

10 

13 

2. 5  7 

10 
13 

2. 75  7 


10 

13 


9co 

(psf) 

124 

105 

81 

66 

79 

159 

149 

156 

214 

185 

152 
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TABLE  11-13  (CONT'D) 


PHASE 

Configu  rat  ion 
C-7a  {VI)(cont'd.  ) 


C-9  (III) 


II  -  PARACHUTE  TEST  SUMMARY 


Max.  Dia.  Mach  qoo 

(In.  )  No.  x/d  (psf) 


6 


8 


2.  9  7 

2.  98  7 

10 

13 

3.  2  7 

10 

13 

3.  5  7 

10 

3. 75  7 

10 
12 

4.  0  7 

10 

4.  3  7 

4.  65  7 

10 

3.  0  5 

7 

10 

3.  2  4 

7 

10 

3.  5  4 

7 

10 

3.  75  5 

7 

10 


134 

125 


105 


81 


66 


158 

149 

157 

123 


105 


81 


66 
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TABLE  11-13  (CONT'D) 


PHASE  11 

-  PARACHUTE 

TEST  SUMMARY 

Max.  Dia. 

Mach 

'loo 

Configuration 

(In.) 

No. 

x/d 

(psf) 

C-3b  (IV) 

8 

3.  0 

10 

123 

3.  2 

6 

105 

9 

10 

3.  5 

7 

81 

10 

3.  75 

7 

74 

10 

C-6a  (VIU) 

8 

2.  5 

7 

185 

8 


2.  75  6  152 

7 

9 


3.0  6  123 

7 

9 

3.2  6  105 

7 

3.  5  6  81 

8 

3.  75  5  66 

8 

C.R.  -2b  (IX)  3. 75  8  79 

10 

12 

10  143 
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TABLE  11-13  (CONT'D) 


PHASE  U  -  PARACHUTE 

TEST  SUMMARY 

Max.  Dia. 

Mach 

^QD 

Configuration  (In.  ) 

No. 

X  /d 

(psf) 

Hemisflo  (X)  8 

2.  3 

8 

214 

10 

1  3 

2.  75 

8 

152 

10 

13 

0 

3.  0 

8 

124 

10 

12 

* 

3.  2 

8 

104 

10 

13 

3.  5 

8 

81 

10 

13 
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contributed  to  the  attainment  of  the  large  amount  of  data  gathered  in 
this  program.  The  greatest  contribution,  however,  was  due  to  the 
favorable  performance  characteristics  of  the  model  canopies,  since 
oscillations  caused  by  improper  design  lead  to  rapid  failure  of  model 
components.  For  example,  one  model  canopy  (C-7a)  survived,  with  no 
apparent  damage,  four  hours  of  operation  in  the  Mach  number  range 
from  M  =  2.  3  to  M  =4.65  in  this  program,  after  having  been  used  in 
the  M  =  2  high  dynamic  pressure  tests  in  the  Phase  la  test  program 
conducted  in  the  Cook  Technological  wind  tunnel  facility. 

Several  other  models  were  subjected  to  several  hours  of  tunnel 
operation  with  either  little  or  no  damage. 

5..  Results  of  Test  Program  and  Data  Analysis 

a.  Parachute  Drag 

As  shown  in  Table  11-13,  three  model  configurations,  name¬ 
ly  C-6a  (I),  C-6a  (V)  and  C-7a  (VI),  were  tested  from  Mach  2.  3 
to  4.  0,  with  the  latter  two  models  tested  to  M  =  4.  65,  Most 
other  models  were  tested  only  to  Mach  3,  75  due  to  limitation  in 
available  tunnel  testing  time.  Drag  coefficients  based  upon 
design  dimensions  are  plotted  versus  Mach  number  in  Figure 
11-28  for  several  configurations,  and  canopy  positions.  In  this 
figure,  it  may  be  noted  that  drag  coefficients  for  most  configura¬ 
tions, reach  a  peak  at  M  =  3.  75  drop  rapidly  after  this  Mach  num¬ 
ber,  and  level  off  at  Mach  numbers  above  4.  0.  A  similarity  to 
these  data  may  be  noted  in  Figure  11-29  for  the  case  of  a  solid 
cone  with  45  degree  half  angle  tested  in  the  wake  of  a  forebody. 
Figure  11-29  is  based  on  data  from  References  4  and  5. 

It  may  be  suspected  in  the  parachute  cases,  that  this  effect 
is  a  reflection  of  either  Reynolds  number  or  dynamic  pressure 
effects.  This  is  strongly  evidenced  by  comparing  the  tunnel 
dynamic  pressures  versus  Mach  number  during  the  conduction 
of  the  Phase  II  program  with  the  associated  drag  coefficients  of 
two  model  configurations,  shown  in  Figure  11-30,  where  maxi¬ 
mum  drag  coefficients  for  model  configurations  C-6a  (V)  and 
C-7a  (VI)  correspond  to  the  minimum  dynamic  pressure  at 
M  =  3.  75.  The  inflated  diameters  of  these  models,  however, 
as  measured  from  the  Schlieren  movies  do  not  appear  to  vary 
significantly  at  M  =  3.  75  from  those  measured  at  other  Mach 
numbers.  For  example,  the  average  inflated  to  constructed 
diameter  ratio  for  test  model  C-7a  (VI)  at  Mach  3.  75  is  0,  89 
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igure  11-28.  Drag  Coefficient  vs.  Mach  Number,  Several  Phase  II  Models 
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and  the  value  of  this  ratio  at  Mach  4.  0  is  0.  90.  A  significant  re¬ 
duction  in  drag  coefficient  would  normally  be  expected  to  be 
associated  with  a  reduction  in  inflated  area.  The  effect  of  Rey¬ 
nolds  number  on  supersonic  parachute  performance  has  not 
generally  been  found  to  be  significant  in  previous  test  programs. 
The  tests  involving  solid  cones  referred  to  above  were  all  con¬ 
ducted  with  Reynolds  numbers  per  foot  in  excess  of  2  x  10^, 
which  would  be  expected  to  represent  turbulent  flow  conditions 
in  these  tests.  Further  study  and  testing  is  required  to  define 
the  nature  of  the  flow  conditions  in  the  Mach  number  range  from 
3  to  4  and  above. 

As  stated  in  Reference  6,  the  qualitative  characteristics 
of  turbulent  boundary  layers  in  supersonic  flows  up  to  Mach  4 
are  similar  to  those  in  incompressible  flow.  At  speeds  above 
this  value,  new  effects  due  to  the  coupling  of  turbulent  velocity 
or  vorticity  field  with  the  fluctuations  in  state  variables  become 
significant.  Since  the  characteristics  of  the  wake  of  a  body  are 
governed  by  the  boundary  layer  properties  of  the  body,  these 
effects  may  correspondingly  become  significant  in  the  perform¬ 
ance  of  a  trailing  decelerator  which  is  submerged,  or  partly 
submerged  in  the  wake  at  high  supersonic  speeds. 

In  considering  the  drag  effectiveness  of  the  tested  model 
canopies,  reference  may  be  made  to  Figure  1-3,  Flow  deviation 
through  the  canopy  shock  w'as  shown  to  result  in  the  spreading  of 
the  wake  after  the  shock  into  a  conical  region  which  extends  to 
the  canopy  inlet.  Within  this  region,  an  essentially  constant 
pressure  region  exists  which  satisfies  pressure  boundary  require¬ 
ments  across  the  normal  portion  of  the  shock  on  the  wake  center- 
line,  and  on  the  diverging  flow  boundary. 

It  may  be  expected  that  the  drag  of  the  canopy  under  these 
conditions  is  equivalent  to  the  drag  on  a  cone  having  the  same 
half  angle  as  that  represented  by  conical  diverged  flow  region. 

In  Figure  11-31  is  shown  the  average  conical  shock  wave  angle  as 
measured  from  Schlieren  photographs  of  Model  C-7a  (VI)  tests. 
Also  shown  is  the  half  angle  cone  which  would  produce  the  equi¬ 
valent  drag  to  that  measured  in  the  tests,  and  the  corresponding 
shock  wave  angles.  Figure  11-31  represents  a  canopy  location  of 
seven  forebody  diameters,  and  Figure  11-32  is  shown  to  repre¬ 
sent  similar  characteristics  at  9.  5  forebody  diameters  aft  of  the 
forebody  base.  Corresponding  measured  and  calculated  shock 
wave  angles  are  shown  for  the  model  configuration  C-6a  (I)  in 
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Figures  11-3  3  and  11-34, 


That  a  stable  flow  condition  may  be  maintained  w'hich 
satisfies  pressure  boundary  requirements  will  now  be  discussed. 
By  the  employment  of  the  methods  to  determine  wake  flow- 
characteristics  discussed  in  Section  I-D,  the  pressure  ratios 
across  the  normal  portion  of  the  canopy  shock  wave  on  the  w-ake 
centerline  are  determined  for  various  x /d  locations  aft  of  the 
forcbody.  These  were  determined  by  the  procedure  described 
in  Section  I  for  various  free  stream  Mach  numbers.  Examples 
of  these  calculations  are  shown  in  Figures  11-35,  11-36,  11-37, 
11-38,  Conical  surface  pressure  coefficients  for  various  cone 
half  angles  at  the  corresponding  free  stream  Mach  numbers  are 
also  shown.  The  agreement  of  the  above  calculated  data  with 
the  performance  of  configuration  C-7a  is  shown  in  Table  11-14 
and  illustrated  in  Figures  11-  35,  11-  36,  11-  37  and  11-38,  InTable 
11-14,  equilibrium  conditions  (examples  of  which  are  shown  in 
Figures  11-35,  11-36,  11-37,  and  11-38.)  are  tabulated  and  corre¬ 
sponding  half  cone  angles  as  determined  by  drag  for  configuration 
C-7a  arc  showTi.  Close  agreement  of  the  calculated  data  is  appar¬ 
ent.  During  the  conduct  of  the  Phase  11  tests,  it  was  found  that 
the  canopy  location  had  a  very  marked  effect  on  the  performance 
of  all  configurations  tested,  and  optimum  canopy  locations 
appeared  generally  to  be  in  region  of  8  calibers  or  less  aft  ofthe 
forebody. 

TABLE  11-14 

CORRELATION  OF  CALCULATED  SHOCK  VERTEX  POSITION  AND 
EQUILIBRIUM  DIVERGED  CONICAL  REGION  WITH  PARACHUTE  DRAG 


M 


2.  75 

3.  0 
3.2 
3.  5 

3.  75 

4.  00 
4.  30 
4.  65 
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MEASUREMENTS 

FOR  MODEL  C 

-7a  (VI) 

Equilibrium 

Equivalent 

Shock  Vertex  Half  Cone  Angle 

Drag  Half  Cone  , 

X  /d 

(Deg.  ) 

(Deg.  ) 

7.  0 

28,  0 

- 

7.  2 

28.  4 

28.  2 

7.  0 

28.  0 

28.  2 

7.  3 

27.  7 

28.  0 

7.  3 

27.2 

29.  3 

7.  1 

26.  5 

27.4 

7.  0 

25.  6 

- 

6.4 

24.  5 

25.  2 

844 

84 

6g  ~Deg. 


Figi. T  11-35.  Wake  Centerline  Normal  Shock  (Vertex)  Pressure  Ratio  Vs,  x/d, 


and  Cone  Surtace  Pressure  Ratio  Vs.  Cone  Half  Angle  at 


Mod  -  3.00 


25  26  27  28  29  Cone  half  angle, 

0  ~  Deg. 

8 

Figure  11-36.  Wake  Centerline  Normal  Shock  (Vertex)  Pressure  Ratio  Vs.  x/d, 
and  Cone  Surface  Pressure  Ratio  Vs.  Cone  Half  Angle  at 
M  00  =  3.  50 
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25  26  27  28  29  30  31  Cone  half 

angle, 

Sg  ~  Deg. 

Figure  11-37.  Wake  Centerline  Normal  Shock  (Vertex)  Pressure  Ratio  Vs. 

x/d,  and  Cone  Surface  Pressure  Ratio  Vs.  Cone  Half  Angle 


25  26  27  28  29  cone  half  angle, 

0  ~  Deg. 

8 

Figure  11-38.  Wake  Centerline  Normal  Shock  (Vertex)  Pressure  Ratio  Vs. 

x/d,  and  Cone  Surface  Pressure  Ratio  Vs.  Cone  Half  Angle 
at  M  00  =  4.  00 
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Upon  examination  of  Figures  11-35  through  11-38,  it  may 
be  noted  that  as  x/d  is  increased,  the  cone  angle  for  equilibrium 
conditions  also  increases.  /Uthough  the  apparent  deflection 
angles  re.Ui/.ed  in  the  examined  cases  are  greater  than  those 
through  an  oblique  shock,  the  high  pressures  (near  stagnation) 
existing  within  the  canopy  cause  the  flow  to  deviate  further.  It 
may  be  anticipated  that  there  will  be  a  limiting  downstream 
canopy  location  for  which  stable  flow  conditions  will  prevail.  This 
is  premised  on  the  assumption  that  the  flow  conditions  described 
earlier  are  responsible  for  the  excellent  canopy  performances 
realized,  particularly  at  low  x/d  canopy  positions. 

A  summary  of  performance  data  for  the  Phase  II  test 
program  is  given  in  Table  11-15.  Drag  coefficients  are  based 
upon  the  actual  constructed  frontal  area  in  all  cases. 

b.  Parachute  Stability 

In  general,  the  models  tested  in  this  program  were  highly 
stable  with  respect  to  the  point  of  suspension.  A  definite  trend, 
however,  was  noted  in  the  stability  of  the  canopy  as  a  function  of 
canopy  location  aft  of  the  forebody.  This  was  discussed  pre¬ 
viously  in  the  subsection  on  parachute  drag.  One  notable  excep¬ 
tion  to  the  generally  stable  performance  characteristics  of  the 
test  models  was  configuration  C-3b  (IV)  which  exhibited  quite 
violent  instability  about  the  point  of  suspension.  Upon  examina¬ 
tion  of  the  model  after  test,  considerable  damage  was  evident  to 
the  roof  mesh  material  at  one  point  near  the  mesh  attachment 
point.  Close  examination  of  the  Schlieren  movies  also  indicated 
signs  of  failure  in  this  area,  but  it  is  impossible  to  determine 
whether  this  damage  occurred  early  in  the  tests,  causing  a 
destabilizing  assymmetry  to  the  vent.  This  model  configuration 
was  a  slightly  modified  version  of  the  Phase  la  Model  C-3a, 
which  exhibited  high  stability  in  the  Phase  la  tests. 

Measurement  of  oscillation  angles  for  the  test  canopies 
was  found  to  be  complicated,  particularly  at  large  x/d  canopy 
positions,  by  the  introduction  of  secondary  oscillations  of  a 
swivel  which  was  used  in  the  tests.  Since  the  weight  of  the  swivel 
(located  at  the  confluence  point)  was  comparable  to  the  model 
canopy  and  suspension  line  weights,  alternate  oscillations  of  the 
two  masses  resulted.  As  a  consequence,  accurate  quantitative 
measurements  of  oscillation  angles  of  the  model  canopies  about 
the  point  of  suspension  were  not  possible,  except  for  the 
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TABLE  11-15 


PERFORMANCE  SUMMARY  OF  THE  PRASE  II  TEST  MODELS 


Mach 

Test  Model 

No. 

x/d 

(averaged^) 

c 

C-6a  (I) 

2.  i 

8.  0 

2.  1 

0.  89 

0.  67 

(d^.  =  7.  9  in.  ) 

9.  7 

0.  5 

0.  87 

0.  70 

2.75 

8.  5 

•1.  2 

0.  87 

0.  64 

9.  5 

1.  5 

0,  82 

0.  63 

12.  5 

5.  0 

0.  85 

0.  57 

3.  0 

8.  0 

6.  5 

0,  67 

0.  56 

9.  5 

7.  5 

0.  81 

0.  54 

3.  2 

8.  0 

10.  3 

0.  91 

0.  61 

9.  5 

7.  5 

0.  87 

0.  55 

3.  5 

8.  0 

3.  0 

0.  79 

0.  63 

9.  5 

8.  0 

0.  85 

0.  57 

3.  75 

8.  0 

7.  4 

0.  82 

0.69 

9.  0 

6.  6 

0.  81 

0.  62 

4.  0 

8.  0 

3.  9 

0.  78 

0.  42 

C-6b  (II) 

2.  75 

8.  5 

8.  1 

0.  84 

0.  54 

(d^^  8.  0  in.  ) 

9.  9 

7.  2 

0.  84 

0.  53 

3.  00 

8.  2 

3.  2 

0.  83 

0.  53 

10.  0 

9.7 

0.  82 

0.  51 

3.  20 

8.  5 

8.  4 

0.  83 

0.  56 

10.  0 

2.  4 

0.  89 

0.  51 

3.  50 

8.  5 

10.  0 

0.  87 

0.  60 

9.6 

5.  0 

0.  86 

0.  55 

3.75 

8.4 

5.  7 

0.  88 

0.  67 

9.  8 

7.  8 

0.  90 

0.  59 

C.R.  -2a  (VII) 

2.  3 

7.  0 

0.  0 

0.  93 

0.  64 

(dc  =  5.40  in.  ) 

10.  0 

4.  8 

0.  93 

0.  79 

11,  8 

2.4 

1.  05 

0.  79 

2.75 

6.  9 

5.  0 

1.  02 

0.  69 

9.  7 

4.  9 

1.  04 

0.  77 

11.  9 

4.  0 

0.  94 

0.  75 
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TABLE  11-15  (CONT'D) 


PERFORMANCE  SUMMARY  OF  THE  PHASE  II  TEST  MODELS 


Test  Model 

Mach 

No. 

x/d 

dm  /d 

(average  d^) 

c 

C.  R.  -2a  (VII) 

3.  00 

6.  7 

2.  9 

0.  94 

0.  67 

(dj,  =  5.  40  in.  ) 

9.  3 

4.  1 

0.  96 

0.  70 

(cont 'd.  ) 

12.  2 

5.  1 

0.  97 

0.  73 

3.  2 

6.  7 

15.  3 

0.  87 

0.  80 

9.  8 

1  0.  0 

0.  92 

0.69 

11.  5 

9.  5 

0.  87 

0.  7  3 

3.  50 

6.  8 

12.  1 

0.  96 

0.  83 

9.  5 

5.  0 

0.  96 

0.  70 

11.  5 

2.  7 

1.  00 

0.  64 

3.75 

6.  6 

15.  5 

1.  00 

0.  89 

9.7 

4.  7 

1.  00 

0.  74 

12.  4 

4.  8 

1.  00 

0.69 

C -6a  (v) 

2.  50 

7.  0 

0.  0 

0.  92 

0.  58 

(d  =  5.  7  in.  ) 

9.  9 

2.  7 

0.  98 

0.  71 

c 

12.  8 

5.  0 

0.  97 

0.  70 

■  7.  0 

0.  0 

0.  88 

0.  58 

Low  , 

9.  3 

2.  3 

0.  91 

0.  73 

q 

,12.  0 

4.  6 

0.  97 

0.  72 

2.  75 

7.  0 

0.  0 

0.  89 

0.  49 

9.  9 

2.  4 

0.  88 

0.  64 

1'2.  5 

3.  8 

0.  93 

0.  64 

3.  00 

7.  0 

5.  1 

0.  96 

0.61 

9.  7 

3.  0 

1.  00 

0.  60 

12.  6 

2.  8 

0.  93 

0.  62 

3.  20 

6.  5 

9.  2 

0.  93 

0.  70 

9.  5 

6.  0 

0.  90 

0.  61 

12.  0 

6.  0 

0.  95 

0.  61 

3.  50 

7.  0 

7.  5 

0.  91 

0.  73 

10.  0 

2.  9 

0.  90 

0.  61 

12.  0 

2.  4 

0.  91 

0.  56 
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TABLE  11-15  (CONT'D) 


PERFORMANCE  SUMMARY  OF  THE  PHASE  II  TEST  MODELS 


Mach 

o 

Test  Model 

No. 

X  /d 

Aq 

(average  d^) 

c 

C-6a  (V) 

3.  75 

7.  0 

2.  7 

0.  93 

0.  85 

(cl^  =  5.  7  in.  ) 

10.  0 

4.  1 

0.  97 

0.  71 

(cont 'd.  ) 

12.  0 

1.4 

0.  97 

0.61 

'1.  0 

7.  0 

2.  4 

0.90 

0.  52 

10.  0 

3.  8 

0.  91 

0.  55 

4.  3 

7.  0 

6.  1 

0.  83 

0.  43 

9.  1 

4.  5 

0.  93 

0.  48 

11.  7 

4.  1 

0.  88 

0.  51 

4.  65 

7.  3 

3.  9 

0.  88 

0.  41 

9.  0 

4.  8 

0.  84 

0.  41 

11.  7 

0.  61 

0.  86 

0.  44 

C-7a  (VI) 

2.  3 

6.  7 

0.  0 

0.  91 

0.  58 

(dj.  =  5.  9  in.  ) 

2.  75 

6.  8 

0.  0 

0.  92 

0.  54 

2.  90 

6.  8 

0.  0 

0.  88 

0.  55 

2.  98 

6.  7 

0.  0 

0.  95 

0.  54 

9.  5 

4.  5 

0.  90 

0.  54 

12.  0 

5.  0 

0.  92 

0.  53 

3.  2 

6.  8 

2.  5 

0.  92 

0.  55 

9.  3 

3.  9 

0.  87 

0.  53 

12.  5 

0.  4 

0.  88 

0.  53 

3.  5 

6.  7 

1.  8 

0.  88 

0.  59 

9.  5 

2.  8 

0.  87 

0.  52 

3.  75 

6.7 

3.  6 

0.  89 

0.64 

9.  5 

1.  8 

0.  90 

0.  56 

12.  0 

6.  0 

0.  92 

0.  53 

4.  0 

6.  8 

0.  0 

0.  90 

0.  52 

9.  5 

3.  8 

0.  90 

0.49 

4.  3 

7.  0 

3’.  4 

0.  80 

0.45 

4.  65 

7.  0 

3.  4 

0.  92 

0.  44 

9.  0 

4.  8 

0.  83 

0.41 
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TABLE  11-15  (CONT'D) 


PERFORMANCE  SUMMAl^Y  OF  THE  PHASE  II  TEST  MODELS 


Test  Models 

Mach 

No. 

X  /d 

[> 

p 

o 

(average  d^^) 

S 

c 

C-9  (III) 

3.  00 

4.  5 

0 

0.  66 

0.21 

(d^  =  8.  0  in.  ) 

6.  6 

0 

0.  68 

0,  29 

9.6 

0 

0.  68 

0.  27 

3.  2  0 

4.  5 

0 

0.  75 

0.  36 

6.  5 

0 

0.  77 

0.  37 

9.  5 

0 

0.  75 

0.  26 

3.  50 

4.  5 

0 

0.  70 

- 

6.  7 

0 

0.  79 

0.  35 

9.6 

0 

0.  77 

0.  27 

3.  75 

4.  5 

0 

0.  67 

- 

6.  8 

0 

0,  74 

0.  38 

10.  0 

0 

0.  78 

0.  32 

C-3b  (IV) 

3.  0 

9.  8 

Model 

Damaged 

0.  32 

(dc  =  7.  7  in,  ) 

3.  2 

5.  8 

Model  Damaged 

0.40 

7.  0 

Model  Damaged 

0.  38 

10.  0 

Model  Damaged 

0.  38 

3.  5 

7.  0 

3.  4 

0.  78 

0.40 

10.  0 

7.  1 

0.  74 

0.  37 

3.  75 

7.0 

3.4 

0.  80 

0.40 

9.  5 

- 

- 

0.  34 

C-6a(VIII) 

2.  5 

7.  0 

11.  8 

0.  68 

0.  39 

(d^  =  7.  9  in.  ) 

7.  5 

1.  3 

0.  67 

0.  37 

2.  75 

7.  0 

0.  7 

0.  63 

0.  24 

7,  5 

8.  5 

0.  56 

0.  25 

9.  8 

6.  1 

0.  71 

0.  30 

3.  0 

6.  9 

0.  0 

0.  71 

0.23 

7.  5 

0.  0 

0.  60 

0.  22 

TO.  1 

1.  7 

0.  68 

0.  24 

3.  2 

6.  7 

5.  3 

0.  69 

0.  27 

3.  5 

6.  1 

0.  0 

0.  81 

0.  32 

7.  9' 

9.  0 

0.  81 

0.  28 
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TABLE  11-15  (CONT'D) 


SPERFORMANCE  SUMMARY  OF  THE  PILASE  II  TEST  MODELS 
Mach 


Test  Models 

No. 

x/d 

t> 

p 

o 

(average  d^) 

C-6a  (VIII) 

3.  75 

4.  7 

0.  0 

0.  69 

0.  18 

(dc  =  7.  9  in.  ) 

6.  6 

1.  8 

0.  69 

0.  28 

(cont'd.  ) 

C.  R.  -2b  (IX) 

3.  75 

8.  0 

7.  0 

0.  88 

0.  69 

(d^  =  8.  0  in.  ) 

10.  0 

4.  8 

0.  83 

0.  54 

12.  0 

4.  0 

0.  81 

0.  52 

Hemisflo  (X) 

2.  3 

8.  5 

0 

0.  43 

0.  10 

(dc  =  8.  0  in. ) 

10.  0 

0 

0.  43 

0.  13 

12.  7 

1.  5 

0.  43 

0.  14 

2.  75 

8.  3 

0 

0.  27 

0.  05 

9.  7 

0 

0.  34 

0.  06 

12.  8 

3.  0 

0.  39 

0.  11 

3.  0 

8.  3 

1.  2 

0.  46 

0.  13 

9.  8 

1.  2 

0.  43 

0.  12 

12.  7 

0.  9 

0.  45 

0.  13 

3.  20 

8.  3 

0 

0.  48 

0.  14 

9.  8 

1.  0 

0.  47 

0.  15 

12.  6 

0 

0.  48 

0.  16 

3.  50 

8.  3 

0 

0.  48 

0.  11 

9.7 

0 

0.  44 

0.  12 

12. .7 

0.  9 

0.  48 

0.  12 
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minimum  x/d  positions  where  the  swivel  location  was  close  to 
the  base  of  the  forebody.  As  pointed  out  previously,  observa¬ 
tions  made  at  the  time  of  the  tests  and  also  from  the  Schlieren 
movies  indicated  the  general  trend  of  increased  stability  with 
decreasing  x/d  canopy  positions  from  an  x/d  of  near  13  to  an 
x/d  of  7.  For  the  six  inch  inflated  diameter  models  with  12  inch 
suspension  lines,  the  minimum  x/d  position  which  could  be 
accomplished  was  slightly  less  than  x/d  r  7.  For  the  8  inch 
inflated  diameter  models  with  16  inch  suspension  lines,  the 
minimum  distance  was  x/d  =  4.  However,  in  this  latter  case, 
configuration  C-6a  (VIII),  violent  instabilities  were  evidenced 
at  x/d  locations  both  greater  and  less  than  the  range  from  x/d  = 

5.  5  to  8.  The  proximity  of  the  large  diameter  canopy  to  the 
wake  trailing  shock  apparently  caused  the  instability  at  x/d  loca¬ 
tions  less  than  5.  5.  Since  the  suspension  line  angle  for  this  con¬ 
figuration  is  approximately  20°  and  the  inlet  cone  angle  is  10°, 
the  suspension  line  inward  load  component  is  apparently  greater 
than  the  net  outward  pressure  differential  at  the  canopy  inlet, 
leading  to  indentations  of  some  of  the  conical  inlet  panels.  This 
was  visually  apparent  during  the  tests  of  this  configuration.  For 
this  short  line  length  configuration,  the  wake  conditions  of  the 
suspension  line  confluence  point  may  also  upset  the  equilibrium 
conditions  described  in  the  previous  subsection.  Large  excursions 
of  the  shock  wave  were  noted  for  this  configuration  with  the  most 
upstream  position  of  the  shock  wave  forward  of  the  confluence 
point.  Interaction  of  the  canopy  shock  wave  with  the  suspension 
line  boundary  layer  could  induce  a  pressure  ratio  across  the 
shock  which  is  sufficient  (Reference  2)  to  separate  a  turbulent 
two  dimensional  boundary  layer,  resulting  in  a  separated  low 
energy  flow  region  extending  to  the  canopy  inlet.  In  Table  11-15, 
are  listed  oscillation  angles  about  the  point  of  suspension  at  the 
low  x/d  canopy  positions,  for  Phase  II  test  configurations. 

c.  Parachute  Inflation 

Inflation  characteristics,  except  for  test  models  C-9  and 
Hemisflo,  were  in  general  good.  In  particular,  the  6  inch  maoci- 
mum  diameter  models  were  outstanding.  In  Figure  11-39  are 
shown  Schlieren  photographs  of  test  model  C-7a  (VI)  at  various 
Mach  numbers,  and  in  Figure  11-40  are  shown  Schlieren  photo¬ 
graphs  of  test  model  C-6a  (V).  In  Table  11-15  are  listed  average 
inflated  diameter  to  constructed  diameter  ratios  for  Phase  11 
configurations.  It  may  be  noted  that  average  inflated  to  construc¬ 
ted  diameter  ratios  of  the  6  inch  diameter  models  are  in  general 
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M  =  3. 75.  x/d  =  7 


M  =  4.  65,  x/d  =  7 


Figure  11-39.  Schlieren  Photographs  of  Phase  II  Model  C-7a(VI)  at  Several 
Mach  Numbers.  Film  Speed:  1,000  fps 
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greater  than  those  for  the  8  inch  models.  Also,  the  percentage 
of  inflation  of  the  8  inch  models  with  16  inch  lines  are  greater 
than  those  for  tlie  8  inch  model  with  8  inch  lines  (C-6a  (VIII)). 

D.  Summary  of  Test  Model  Performance 

Those  configurations  which  were  characterized  by  a  high  ratio  of  coni¬ 
cal  inlet  diameter  to  maxirnuir.  diameter,  a  low  cone  inlet  angle,  a  flat  roof, 
with  low  inlet  cone  porosities  and  relatively  high  roof  porosities  were  out¬ 
standing  in  performance  characteristics  in  the  Phase  la  and  Phase  II  wind 
tunnel  test  programs.  These  configurations  (such  as  models  C-6a,  C-6b, 
C-7a,  C-7b,  CR-2a  and  CF-Zb)  also  represent  the  lowest  ratio  of  total  cloth 
area  to  projected  frontal  area  (2.  0).  of  all  models  tested.  This  ratio  is  also 
lower  than  conventional  circular  flat  canopies  which  normally  have  a  ratio  of 
2.25  in  the  speed  ranges  where  they  perform  satisfactorily. 

In  the  design  of  the  models  listed  above,  a  wide  range  of  flat  roof  con¬ 
structions  was  employed.  Total  porosity  ranges  of  these  models  were  from 
8  to  18  percent.  Corresponding  roof  porosities  ranged  from  14  to  34  per¬ 
cent,  and  inlet  porosities  for  all  of  the  aforementioned  models  were  low 
(2  percent). 

Very  little  difference  was  noted  in  the  relative  performance  of  these 
configurations,  although  in  general,  those  which  exhibited  the  highest  stability 
about  the  point  of  suspension  were  constructed  with  the  highest  total  porosity 
(e.  g. ,  C-7a).  This  appears  to  confirm  early  data  with  respect  to  the  effect 
of  total  porosity  on  stability.  As  pointed  out  previously,  however,  these 
differences  were  small  even  though  a  relatively  wide  range  of  total  porosities 
was  employed  in  the  model  designs. 

The  configurations  C-6,  C-7,  and  CR-2  resulted  from  a  series  of 
modifications  of  the  Phase  I  test  models  A  and  B. 

Although  the  model  C-9,  which  evolved  from  the  Phase  I  model  D,  did 
not  perform  outstandingly  in  the  Phase  la  and  Phase  II  test  programs,  it  is 
felt  that  a  considerable  improvement  in  the  inflation  characteristics  of  this 
configuration  would  result  if  the  ratio  of  inlet  diameter  to  maximum  diame¬ 
ter  were  increased.  This  particular  configuration  in  such  a  modified  version 
would  be  particularly  adaptable  to  those  applications  where  extreme  stability 
of  the  trailing  device  is  essential.  Since  the  drag  of  a  parachute  is  proportion¬ 
al  to  the  inlet  area,  a  compromise  of  required  drag  and  stability  could  be 
attained  with  this  configuration  by  an  appropriate  selection  of  the  ratio  of 
inlet  diameter  to  projected  diameter. 
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The  Hemisflo  c onfif»uralion  which  was  tested  in  the  Phase  I  and  Phase  11 
test  programs  failed  to  inflate  properly  at  any  time  in  the  speed  range  from 
M  =  2.  3  to  3.  75.  Since  the  high  inflation  characteristics  of  the  Models  C-6, 
C-7  and  CR-2  were  primarily  attributed  to  the  low  porosities  employed  in  the 
inlet  region,  as  well  as  the  low  total  porosities  used,  it  is  not  surprising 
that  this  configuration  did  not  inflate  well.  This  is  particularly  striking  since 
hemisflo  parachutes,  like  most  conventional  ribbon  constructions,  are 
designed  with  the  maximum  localized  porosities  in  the  inlet  and  skirt  regions 
and  are  generally  constructed  with  relatively  high  total  porosity.  It  is  not 
expected  that  a  mere  reduction  in  total  porosity  of  this  canopy  would  result 
in  its  satisfactory  performance,  since  a  porosity  distribution  sucii  us  to  pro¬ 
vide  low  values  in  the  inlet  area  have  been  demonstrated  in  this  program  to 
be  desirable. 
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SECTION  III 


SYMBOLS  AND  NOTATION 

A  Wakf  pcTranictc  r 

a  Speed  of  sound 

a,  b,  c  Constants  for  streamline  fitting  in  Section  III-E 
B  Wake  parameter 

CpjCy  Specific  heals  at  constant  pressure  and  constant  volume, 
respectively 

d  Diameter  of  body 

f  Function  representing  entire  streamline 

K  Correction  factor  for  boundary  conditions  at  shock  wave 

k,  1,  m  Constants  for  normal  fitting  in  Section  III,  E 
Reference  length 
M  Mach  number 

n  Arc  length  of  normal 

P  Point  in  flow  field 

p  Pressure 

Q  Infinitesimal  quantity 

R  Gas  constant 

r  Number  of  points  used  in  least  squares  fit 

S  Entropy 

s  Arc  length  on  streamline 
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SECTION  III 


SYMBOLS  AND  NOTATION 
(CONT'D) 

T  Temperature 

V  Velocity 

X,  y  Cartesian  coordinates 

a, 13  Orthogonal  curvlinear  coordinates 

y  Ratio  of  specific  heats 

S  Deflection  angle  of  streamline  through  shock  wave 

€  Constant 

^  Infinitesimal  quantity 

T)  Function  representing  entire  normal 

0  Angle  of  flow  direction 

p  Density 

St  ream  function 
w  Shock  wave  angle 

SUBSCRIPTS 

i  Streamline  index 

j  Normal  index 

s  Quantities  immediately  behind  the  shock  wave 

o  Total  flow  conditions  of  undisturbed  flow 

00  Freestream  conditions 

1,2  Quantities  before  and  behind  shock  wave,  respectively 

SUPERSCRIPTS 
h  Level  of  iteration 
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SECTION  III 


ANALYTICAL  STUDY 

The  purpose  oi  the  theoretical  phase  of  this  program  was  the  establish¬ 
ment  of  a  procedure  whereby  feasible  parachute  shapes  could  be  predicted 
analytically. 

In  order  'n  accomplish  the  above,  it  w^as  necessary  to  decide  upon  a 
technique  or  techniques  for  solving  the  flow  equation  for  all  regions  down¬ 
stream  of  the  bow  siiock  wave.  One  approach  would  have  treated  each  of 
the  subsonic,  transonic,  and  supersonic  regions  individually,  thereby  re¬ 
quiring  three  diflerent  analytical  techniques.  A  second  approach  was  one 
whereby  one  technique  would  treat  all  three  regions.  In  addition,  this 
approach  would  correct  the  shock  shape  and  shock  standoff  distance.  Because 
of  the  provision  rendered  for  shock  adjustment  and  the  obvious  advantage  of 
one  rather  than  three  different  analytical  techniques,  the  second  approach 
was  adopted. 

The  approach  adopted  is  the  method  of  flux  analysis  due  to  Uchida  and 
Yasuhara.  This  is  a  method  for  finding  exact  numerical  solutions  for  the 
direct  problem  in  which  the  shape  of  a  blunt  body  in  a  uniform  supersonic 
flow  is  given  and  the  shape  and  location  of  the  detached  bow  shock  and  field 
behind  it  are  unknown. 

The  method  was  modified  appreciably  so  that  it  would  provide  the  initial 
approximation  to  the  flow  field  as  well  as  successive  corrections  to  it.  In  the 
modification,  some  features  of  the  inverse  problem  were  utilized.  The 
inverse  problem  is  one  whereby  the  shape  of  the  shock  wave  is  given  and  the 
shape  of  the  body  and  details  of  the  flow  are  unknown. 

The  approach  was  entirely  mechanized,  however,  it  was  found  to  pre¬ 
sent  certain  difficulties  impractical  to  circumvent  within  the  time  of  this 
program.  No  reasonable  way  was  found  to  make  compatible  elements  of  the 
inverse  and  direct  approaches. 

Consequently,  recourse  was  taken  to  a  semigraphical  method  parallel¬ 
ing  that  indicated  in  Reference  7,  which  is  described  below. 

A.  Shock  and  Compressible  Flow  Relations 


1.  Uniform  Upstream  Field 

For  the  case  of  the  uniform  upstream  flow  the  shock  conditions 
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necessary  in  the  field  and  shock  calculations  are  obtained  from  the 
standard  shock  and  flow  relations  as  follows: 


For  a  given  environment  and  vehicle  Mach  number  the  following 
parameters  are  known: 

^1‘  '^l*  ^1 

Using  these,  the  remaining  parameters  can  be  calculated: 


Tq  =  Tj  (1  +  o.zMj^)  m-1 

Vj/ao  =  (1  +  0.2Mj2)‘^/^  III-2 

P2/p  1  =  (TMi^sin^oj-  l)/6  III-3 

PZ  =  Pi  'P^/P  i  ) 

Pi  I  Pi  =  6Mj^  sin^w/iMj^sin^u  +  5)  in-5 

Pi  ^  P\  {  Pi  IPi  )  III-6 


M- 


36M2'^sin^w  -  5(M^  ^sin^u  -  1  )(7  M|^  ^sin^u  +  5) 


/  in-7 


(7Mj^^  sin  -  1 )  (Mj'sin^w  +  5) 


T^/Tj  =(7Mi^sin^a)- l)(Mi^sin^u+  5)/36M^^3in'^w  HI-S 


2„,-  2 


2„._2. 


V,/ 


I'  ^o 


M 


2  {T2/Tj)(Ti/tj 


1/2 


III -9 
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1  /tan  S  =  tan  w 
2.  Non-Uniform  Upstream  Flow. 

In  the  case  of  the  analysis  of  the  flow  field  around  a  parachute, 
a  portion  of  the  upstream  flo\^  will  necessarily  be  non-uniform  due  to 
the  presence  of  the  vehicle  wake.  Therefore,  the  analysis  of  the  pre¬ 
vious  section  must  be  n'.odified  and  extended  in  order  to  accommodate 
this  type  of  upstream  field. 

The  equation  for  the  velocity  distribution  in  the  wake  at  a  given 
caliber  distance  aft  of  the  vehicle,  as  discussed  in  Section  I,  is  re¬ 
written  for  convenience  in  this  analysis  and  expressed  as  follows: 


l.ZMj  /(Mj  sin  w-  1) 


III-IO 


V^/ao  =  Vco  /ao  I  1  -  A 


B  - 


where 


Ul-11 


and  A  and  B  are  wake  parameters 

The  resulting  Mach  number  distribution  is: 


Ml 


(Vl/ao) 


1-0,  2  (Vj/a 


■1/2 


III- 12 


The  assumption  is  now  made  that  the  static  pressure  in  the  wake 
after  the  wake  throat  is  constant  and  equal  to  the  free  stream  static 
pressure. 

With  pressure  and  velocity  distributions  known,  it  remains  only 
to  determine  the  distribution  of  density.  Since  this  method  of  flux 
analysis  requires  the  field  have  constant  total  energy,  the  density  dis¬ 
tribution  will  be  calculated  using  this  criterion.  Physically,  this 
criterion  is  applicable  since  the  presence  of  a  vehicle  and  vehicle  bow 
shock  do  not  appreciably  alter  the  total  constant  energy  possessed  by 
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the  uniform  Clow  upstream  of  the  vehicle  bow  shock.  For  a  perfect 
gas,  this  condition  of  total  constant  energy  is  satisfied  by  requiring 
tne  total  temperature,  T^,  to  be  constant. 

Thus: 


Po  =  Po  /RTq 


III-l  3 


whe  re 


Po  =  Pi  (1  +  0.2 


IU-14 


Finally 


III-15 


The  shock  conditions  are  now  obtained  from  equation  III-3 
through  III-l  0. 

B.  Field  and  Shock  Calculations 


1.  Field  Calculations 

The  over-all  scheme  of  this  analysis  is  similar  to  many  numeri¬ 
cal  techniques,  that  is,  to  set  down  an  initial  approximate  solution  to 
the  downstream  field  and  shock  and  then  correct  both  by  successive 
iterations. 

With  the  field  nomenclature  as  shown  in  Figure  III-l,  the  pro¬ 
cedure  for  correcting  the  initial  field  is  as  follows.  Isolate  the 
j  and  j  +  1  normals:  (See  Figure  UI-2). 

On  the  j  normal  boundary  conditions  to  be  satisfied  are 


l//G  =  0  at^=  0  III-l  6 
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(y)/2)  ^  >'1 


at  /S  =  9 


III-17 


or  in  non-dimensional  form 


'Ps 

JJ^ITTiT) 


'-1 

Po^-i 


(y,  /  ■< 


UI-18 


where  €  -  0  in  two  dimensional  flow  and 

e  =  1  in  axially  symmetric  flow- 

The  condition  of  momentum  conservation  through  the  shock  is 
temporarily  relaxed  although  all  conditions  arc  to  be  satisfied  in  the 
final  solution. 


For  a  given  fi.  ld  pattern  and  distribution  of  entropy,  it  is  possi¬ 
ble  to  determine  the  velocity  distribution  along  a  normal  by  the  follow¬ 
ing  equation: 


^i  +1^^^ 


1 


X-1 


III- 19 


where: 


As  =  distance  between  two  successive  normals  (j  and  j+  1) 
measured  eilong  a  streamline 


(Sj  -  S.  ^  ^  )  =  entropy  change  between  the  i  and  i  +  1  stream¬ 
lines 
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The  change  in  entropy  between  streamlinett  is  calculated  from 
the  assumed  shape  of  the  shock  in  the  following  manner: 


■  ^i  +  1 )  /  ^p 


III-20 


Knowing  the  velocity  distribution,  one  can  calculate  the 
mass  flow  distribution  by  tlie  following  equation: 


-(Si  -Sj).  /R 

a  .  . 
o  1,  j 


Xll 

1 


Z  l/y-l 
I  m-21 


whe  re: 


m-22 


Equations  III-19  and  III-21  are  utilized  in  an  iterative  procedure 
in  order  to  satisfy  the  boundary  condition. 


When  the  upstream  flo'y-1-a^-uviir'.^i-ffrythe  right  hand  side 
of  equation Til- 2 3  can  be  integrated  into: 
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The  steps  in  detail  are: 


(1 )  Assume  a  velocity  at  the  shock  wave  j 


Use  equation  III-19  to  compute  V.  .  the  velocity  at  the  point 

1  ^  ^ 

(x,y)i^j.  Knowing  V.  j,  again  use  equation  III- 1 9  to  compute 
the  velocity  V.  j  j,  at  the  point  (x,  y)-_j  •.  Continue  along  the 
j  normal  until  al*l  velocities  are  computed. 

Use  equation  III-21  to  compute  the  mass  flow  at  each  point  along 
the  j  normal 

Knowing  the  mass  flow  along  the  j  normal,  numerically  integrate 
equation  III-23 


"s// 

f  JJL 


)  ^+1  Vo 


the  boundary  condition  is  satisfied.  If  not,  vary  V  and  repeat 
the  above  from  Step  2  until  the  above  criterion  is  satisfied. 

The  distribution  of  the  stream  function  along  the  j  normal  at  the 
given  ^  lines  (streamlines)  is  found  by 


^Vo 

2 


nj// 

f  II 

J  Po\ 


and  the  following  type  relation  is  obtained. 


Po’^o  i(JJ 
2 


III-25 


Figure  III- 3.  Flux  Distribution  Along  A  Normal 
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The  location  of  the  streamlines  is  corrected  by  entering 
Figure  111-3  with  the  stream  function  vidues  of  the  streamline  and 
finding  the  corrected  n/^  location  along  the  normal. 

This  procedure  is  followed  on  every  normal  and  the  new  stream 
lines  are  obtained  by  connecting  the  newly  obtained  corresponding 
points  on  adjacent  normals.  New  normals  are  drawn  thus  complet¬ 
ing  the  present  field  iteration.  This  new  Qow  pattern  will  serve  as 
the  co-ordinate  axes  lor  the  next  iteration. 

For  a  given  shuck  shape  iJie  above  iterative  procedure  is  re¬ 
peated  until  a  consistent  field  solution  is  obtained. 

2.  Shock  Calculations 

Fulfillment  of  the  boundary  condition  requires  that  the  shock 
stream  deflection  angles  coincide  with  the  streamline  angle  immedi¬ 
ately  behind  the  shock  obtained  from  the  field  calculations.  That  is: 


This  condition  will  not  be  satisfied  at  every  shock  point  unless  a  solu 
tion  has  been  attained. 

A  parameter  indicating  the  completeness  of  the  solution  is 
defined  as; 


K 


max 


/8 


max 


m-26 


This  parameter  is  used  to  calculate  new  values  of 
next  shock  iteration  as  follows; 


8 


for  the 


8=  %/K 


in-27 


When  K  =  1  is  obtained,  the  shock  condition  is  completely  satisfied. 
This  condition  will  only  be  attained  when  the  correct  shock  stand-off 
distance  has  been  found. 
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The  dfscription  that  follows  is  the  iterative  scheme  for  shock 
adjustment  that  is  employed  for  a  given  shock  stand-off  distance.  The 
values  of  K  obtained  from  successive  iterations  will  approach  a  limit. 

If  this  limit  value  is  not  K  =  1,  the  shock  stand-off  distance  is  varied 
and  the  entire  iterative  scheme  of  shock  and  field  calculation  is  repeat¬ 
ed.  The  new  stand-off  distance  is  obtained  from  the  previously  ob¬ 
tained  value.  If  K  (in  the  limit)  <  1  the  stand-off  distance  is  decreased. 
If  K  >  1  the  stand-off  distance  is  increased.  This  procedure  is  pur¬ 
sued  until  a  value  of  K  =1  is  obtained.  This  will  be  achieved  only  when 
the  correct  shock  stand-off  distance  is  located. 

The  steps  for  shock  adjustment  in  detail  are: 

(1)  Upon  completion  of  a  field  iteration  the  values  of  w.  and 
{0s)j  are  tabulated. 

(2)  From  equation  UI-lO  obtain  the  values  of  §.  corresponding 
to  the  w,  of  step  (1). 

(3)  Sur%ey  the  lists  of  S  and  (  6  ).  locating  the  maximum  of 

®  i  /ft  i. 

each  and  form  the  ratio  K  =  ^  s'i  max. 

‘^i  max. 

(4)  Obtain  a  new  value  of  B  from  equation  III-27. 

(5)  For  each  value  of  S  of  step  (4)  use  equation  III-IO  to  ob¬ 
tain  the  corresponding  values  of  w.. 

(6)  Draw  a  new  shock  wave  starting  the  axis  at  the  chosen 
stand-off  distance.  Each  segment  of  the  shock  wave  is 
drawn  at  the  appropriate  angle  starting  at  the  inter¬ 
section  of  the  shock  wave  segment,  drawn  from  i  -  1 
streamline,  and  the  i  th  streamline. 

(7)  The  shock  wave  now  intersects  the  streamline  at  different 

locations.  Therefore,  measure  the  new  values  of  (  S  ). 

'  s'l 

and  repeat  procedure  from  Step  3. 

This  procedure  is  repeated  until  K  approaches  a  limit. 
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c. 


Circular  Cyliiulor  ProbK-ms 


To  obtain  inaij^lu  to  ilu-  different  aspects  of  the  problem,  the  calcula¬ 
tion  presented  in  Reference  (7)  "A  Circular  Cylinder  in  Two  Dimensional 
Flow",  was  performed. 

The  correct  shock  stand-uff  distance  was  obtained  from  Reference  (7). 
For  this  stand-off  distance,  as  shown  in  Fij^ure  III--1,  the  third  shock  and 
field  iteration  compared  favorably  with  the  fourth  iteration  as  given  in 
Reference  (7). 

Because  of  this  agreement  and  also  because  the  final  solution  as  ob¬ 
tained  in  Reference  (7)  did  not  differ  appreciably  from  the  -Ith  iteration,  the 
calculation  was  terminated. 

This  calculation  was  done  as  a  semig raphical  procedure.  That  is,  the 
field  calculation  was  performed  on  a  high  speed  digital  computer  up  to  the 
point  where  the  location  of  the  corrected  streamlines  had  been  obtained. 

The  drawing  of  the  new  streamline  and  normals  was  done  graphically 
as  was  the  adjustment  of  the  shock  wave. 

D.  Parachute  Problem 

A  limited  effort  was  expended  in  order  to  obtain  the  flow  field  for  a 
three-dimensional  parachute  in  non-linear  flow;  the  non  linearity  resulting 
from  the  presence  of  the  vehicle  wake  upstream  of  the  parachute. 

In  this  case  only  the  field  downstream  of  the  shock  and  upstream  of  the 
parachute  was  considered.  Hence,  the  shock  wave  obtained  from  Schlieren 
photographs  was  assumed  fixed.  (In  order  to  adjust  the  shock  shape,  the  en¬ 
tire  field  must  be  considered.  ) 

Field  calculations  were  stalled  when  difficulties  were  encountered  in 
that  the  governing  flow  equations  were  not  able  to  successfully  correct  por¬ 
tions  of  the  initial  downstream  field.  This  difficulty  was  not  encountered  in 
the  two-dimensional  cylinder  in  uniform  flow. 

Further  study  would  be  required  to  determine  cause  of  the  difficulties. 

E.  Complete  Automation 

In  order  to  improve  on  the  accuracy  of  the  field  calculation  and  to  re¬ 
duce  the  time  involved,  a  method  has  been  evolved  for  automating  the  entire 
field  calculation.  It  is  presented  to  indicate  the  feasibility  of  complete  auto¬ 
mation. 
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1. 


Nonitncl.il  nr  t 


‘I’lit  followiiifi  norncnclatnrt  will  be  used  in  tin*  construclion  of 
St rtanilint  s  and  norni.ils.  (Stt  Figure  III- 5) 

Let  fj  rtprtstnl  the  tnlirt  i^^'  slrtamlint  and  represent  the 
entire  norni.il . 

Pj,j  shall  be  the  corrected  points  used  to  construct  the  new 
St  reamline . 

P  sh.ill  be  points  on  the  new  streamline, 
n 

P  .  .  shall  denote  the  intersection  of  the  new  streamlines  with 

1'  J 

the  new  normals. 


The  points  P^,  j  lie  on  the  axis  or  on  the  body.  The  points 
Pjij,  where  i  =  j,  lie  on  the  shock  wave,  i.e.,  P.,  .. 

Let  f.  ,  represent  the  second  order  curve  passing  through  the 
points  Pi.l'-'Pi,  j,,.  Pj, 


Let  f-  •  represent  the  second  order  curve  passing  through  the 

I  L  J 

point  P  .  .  (j  >  i  +  2)  and  satisfying  a  least  squares  fit  through  the 

points  Pj,  j  +  1'  Pi,  j  +  2» . .  ^i,j  +  r  ^  num¬ 
ber  of  points  used  in  the  least  squares  fit.  It  will  represent  the 

streamline  in  the  neighborhood  of  the  point  from  which  it  starts, 

^  i.  j’ 


Streamline  Construction 


As  described  previously,  new  streamlines  are  constructed  by 
drawing  a  curve  through  corrected  points.  Because  of  the  scatter 
that  occurred  in  the  corrected  points  for  most  portions  of  the  field, 
passing  a  curve  exactly  through  the  points  would  result  in  an  extreme¬ 
ly  rough  curve.  Therefore,  a  least  squares  curve  technique  will  be 
used  to  fair  in  these  streamlines. 


The  only  region  where  streamline  fitting  through  the  exact 
points  will  be  utilized  is  in  the  neighborhood  of  the  shock  wave  where 
the  correction  to  the  streamlines  are  fairly  uniform. 


» 


> 
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A  quadratic  function  is  used  for  the  streamline  fit  since  the 
fitting  of  the  normals  would  become  extremely  complex  if  higher 
orders  were  used  (with  the  quadratic  it  is  necessary  to  solve  four 
non-linear  simultaneous  equations). 

Thus,  to  start  the  procedure  of  fitting  the  streamlines,  a  quadratic 


fi.i  :  y  =  x-  *  ;  X 


X  +  c. 


1. 1 


111-28 


is  passed  through  three  successive  points. 


P  P  P 

i,  i’  i,  i  +  1*  i,  i  +  f 


The  coefficients  of  equation  111-28  will  be  determined  as  follows: 


y..  =a  X  +b.  .X  +c 

t  1,1  i,  i  1, 1  i,  i  i,  i 


y..,,=a.  .X  ;i.i  +  t)..x..  ,tc,. 

^  1, 1  +  1  1,1  1, 1  +  1  1, 1  1,1  +  1  1,1 


+  b  .  .  X  .  .  ,  +  c 


m-29 


y.  .,-,=a,,x  ..  _+b.  .X,,  ^  +  c,. 

1, 1  +  2  1, 1  1,1  +  2  1, 1  1, 1  +  2  1, 1 


Thus 


1. 1 


X  .  .  X  .  , 

1.  1  1. 1 


^i,i  +  l  i  +  1  ^ 


2 

X  . 


i,i  +  2  i,  i  +  2 


1II-29.1 
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y  .  .  X  •  ■ 

i,  1  1. 1 


y.  .  ,  X  1 

W  ‘  +  1  i, :  +  1 


'i,  i  +  2  i,  i  +  2  * 


III-29.  2 


X  .  .  y  .  .  i 

1,  I  ^  i,  I 


•  .  v 

».  I  +  1  i  +  1 


b  .  . 
!•  1 


i,  i  +  2  y  i,  i  +  2 


III-29.  3 


X  .  . 

1.  1 


X  .  .  y  .  . 

If  1  ^  1,  1 


i,i  +  l  ""1,1  +  1  y  1,1+1 


X  X  y 

1, 1  +  2  1, 1  +  2  1, 1  +  2 


m-29.4 


A  similar  curve  will  be  obtained  for  ^  the  next 

lower  streamline  ' 


f  :y  =  a  x+b.  ,.,x  +  c.  , 

i-1,  i-1  i-1,1-1  1-1, 1-1  1-1,  i-1 


III- 30 


From  the  point  P.  •  ,  ,  working  downstream,  the  stream- 

1  T  w 

lines  will  be  fitted  by  a  least  squares  technique. 
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Ill  order  to  maintain  some  relationship  between  the  P.  .  and 
P*  points  let  x'  ;  i  =  X  j  1. 

i,j  ‘‘J  ''J  , 

Two  conditions  are  to  be  satisfied  at  each  P  ^  .  point;  con¬ 
tinuity  and  continuous  first  derivative.  ’ 

The  problem  is  to  determine  f ■  j . 


X  .  .  =  X  . 

i.J  >.J 


m-31 


I 

y  .  •  =  f  .  .  ,  (-X  .  •  Y  . 

M  i.j  -  1  ^  i,y 


m-32 


-  1  i.j  J  y  i,  ^ 


m-33 


F  .  . 


i.j  +  1  ‘“y  [’‘i.j  tl  ■  '"i.  j  i,j  +  1,  ’'i,j  +  l>] 


It  is  now  required  to  determine  f  •  :  which  is  of  the  form 

J 


*  ^ 
f  .  .  (x)  =  a.  .  X  +  b  X  +  c 
1.  J  i,j  i.J 


m-34 


for  which 


f*.  .  .)  =  £*.  .  ,  (x,  i  ) 

1,J  I.J  1,J-1 


in-35 


^  .  .(x.  .)  =  f  i  i  1  (x  .  .) 

1.  J '  1.  J '  1.  J  -  i  '  i,j' 


m-36 


and  the  sum  of  squares  of  the  deviations 
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■^i.J  = 


_  y  , 

j  ^  3=1  i.  j  j  +  3^  "  y  i,  j  +  3  ] 


is  minimized 


From  III- 3-1  and  UI-J5 


III-37 


1  ;  :  t  i  )  =  a;  i  x^;  ;  +  b.  •  X.  .  +  C-  ■ 

J  ■  1  W  J  1|  j  >»  J  1,  j  1,  j  1,  J 


m-38 


Taking  the  derivative  of  III-34  and  using  III-36  renders 


^  i  i  1  {5<  ;  ;)  =  2  a.  .  X  .  .  +  b.  . 


m-39 


From  III-39 


b.  .  =  f  (x  .  .  )  -  2  a,  .  X 

J  J  -1  i.J  1.  J  1,  j 


III- 40 


Substituting  III-40  into  III-38  renders 


•  L 

^  i  i_i •  )  =  C.  .+  X.  .  f  ,  .  (x.  .)_a.  .  X  .  . 

^  i.J  i.J  i.J-1  ^,j  i,J 


Solving  for  c 


^  ^  .  .  Ax.  .)-x  .  .  f  .  .  ,  (x.  .)+  a.  .  X  .  . 

i.j  i.j-1  hr  i.j  i.j-r  i,j  i,j 


in-41 
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Substituting  III-40  and  into  111-3*1  gives 


*  '  i.j-1  '  i.j-1 


m-42 


From  III- 37  since  a  .  .is  the  only  variable  of  the  function  to  be 
^  •  I  ‘'J 

dete  mimed 


IU-43 


From  III-42 


h.jKj+al]  ==‘\jts  -  + 


m-44 


Substituting  III-42  and  III-44  into  III-43  gives 


I 

Z/  ,  (x.  .  )-  X.  .  f  .  .  ,  (x.  .  )  +  a-  •  x^*  •  + 

S  =  1  \  1|  j  1,  J  ly  J  J  i»  J 

^i»j+s  j  ^i,  j  ]  ^  ii  j  +  s  ”^i,j+s  } 

r  2  2 

tx..-2x  X  +x  r=0 

’■  J  i,  j  ii*  j  +9  i»  j  +  s 
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Solving  for  a-  j 


‘»j  '  +  'i,j  +  s  ’  i>j  ^  i.  j-1  '  i,  j  ]' 

s  =  1  ^ 

HI- 

.  -2x  .  ,  X  .  +x  .  ll  / 

L  ‘.J  i,J  i.J  +  3  ‘.J+sJJ/ 


A  r  2 


y  [x  .  -  2  X  .  .X  .  .  +x^  1 

•;)  »,J  l,J  +  3  l.J  +  sj 


3  =  1 


After  a  •  is  found, 
J 

III-41. 


b  .  and  c  .  are  found  from  III-40  and 


3.  Normal  Construction 


Usually,  the  construction  of  normals  will  begin  at  the  shock; 
however,  a  general  field  point  will  be  discussed  since  the  boundary 
point  differs  only  slightly  procedure. 


The  equation  for  the  normal  segments,  '’7.  will  be  expressed 
as  X  =  f  (y)  since  at  the  axis,  the  normality  condition  would  require 
that  dy/dx  —  oo  while  in  the  form  suggested  dy/dx  =  0 


Thus 


X 


+  •  y  +  m 


ni-46 


To  determine  the  coefficients,  the  following  three  conditions  will 


be  used. 

(1) 

intersects 

f  .  .  at 

'i|  J 

hi 

(2) 

’'i.j 

,  fV  , 

^11 

at  P  .  . 

1.  J 

(3) 

II 

’’i.j 

j_  at  some 

point  P  . 

1- 

to  be  determined. 
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The  first  condition,  V  \  : 
II  h  •  J 

stituting  X  j  j  and  y  j  j  for 


intersects  f  .  . 

i.  J 

X  and  y  in  III -do 


is  obtained  by  sub¬ 
yielding. 


II 


X 


i.  j 


'7 

y  .  . 

J 


+ 


y  .  .  +  m.  . 
i,  j  J  J 


III-47 


The  second  condition,  17.  .  X  f  .  .  at  P".  .  is  obtained 

‘.J  WJ  WJ 

from  111-34  and  III-46  by  determining  the  slopes  of  the  two 

curves  at  the  point  P".  .. 

*  J 


F rom  III  -  34 


From  III-46 


'•J 

=  Za  .  X  +  b.  . 

dx 

dx 

J  .1 

dx 

--^k.  ,y+/.  . 

dy 

dy 

J  J 

The  normality  condition  is  satisfied  by  setting 


<< 


dx 


dy 


m-48 


and  substituting  x' 


and 


for  X  and  y. 


or 


Za.  .  x”,  .  +  b.  .  =  -  (2  k;  ;  y"i  ;  +  ■  ■) 

1,  J  1,  J  J  W  J  '  J  1,  J ' 


in-49 


The  third  condition  j  would  be  met  by  a  condition 

similar  to  111-49.  However,  the  point  of  intersection  is  not  known 
and  this  becomes  a  fourth  variable.  A  fourth  condition  will  be  ob¬ 
tained  by  the  generalized  condition 
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4) 


Ct 

intersects  f  .  ,  . 


Substituting  the  expression  for  ^  j  j  III-46  gives 

2  2 
X  =  k.  .  (a;  1  ;  X  +  bj  ,  ;  X  +  c.  ,  .) 

J  l-l.  J  J  1-1,  J 


III-50 


>1' 

Differentiating  the  expression  for  f 
the  condition  of  11-48  renders 


and  III-46  and  satisfying 


111-51 


Substituting  f  into  III-51  gives 

-  1.  j 


2  a.  .  .X  +  b  1  ;  =  -  2k.  .(a.  ,  .  x  +  b,  .x  +  c.  •)■/•  • 
i-i»J  1  -i|J  i,J  i-l,J  i-l,j  i-l,J 


Ul-52 


Equations  111-47,  111-49,  111-50  and  111-52  are  four 
equations  in  the  four  variables,  k.  .1.  .  rn.  .  and  x  where  x  is 

i.j,  1,  J,  i,  J 

actually  x".  ,  .  one  co-ordinate  of  the  intersection  of  77:  : 

,  ,«  1-1,  J  It  ,  .  ,  ,  >  r,,,  •  J  /  *  J. 


and  f 
P' 


i-i.r 

which  is 


i-l.j 

use 


is  obtained  from  f 


This  determines 


4^  ” J 
for  extending  the  normal  to  the  next  stream- 


Because  of  the  non-linearity  of  the  system  of  four  equations, 
simultaneous  solution  would  be  intractable  .  Instead,  an  iterative 
technique  will  be  applied. 

Unfortunately,  there  is  no  method  of  attack  that  will  guarantee 
convergence.  The  best  approach  is  to  select  initial  values  that  will 
be  close  to  the  solution.  In  the  system  to  be  solved  it  is  fortunate 
that  not  all  equations  contain  all  variables  and  it  will  be  possible  to 
start  the  iteration  by  estimating  two  of  the  parameters. 


ASD-TDR-62-844 


123 


The  parameters  tliat  can  best  be  estimated  are  the  second 
derivative  of  Tj  ,  i.i-.,  k  ,  intersection  co-ordinate 

'  i ,  j '  •  '  ! ,  j 

For  these  estimates  we  use  P;  ,  i  ;;  P  •  and  P.  ,  .  where 

‘  J  i»  J  i-l,  J. 

the  superscripts  indicate  the  iteration  order. 


=  (k.  j'" 

^  V  i+I.j'^i.j  ■  yi  -l,j 


(k.  .)  = 

1.  J 


i+1.  j  ■  ’^i,  j 


Wj  i  -l>j 


^^i+l,j  ^i-ii  jM  ^i+1,  j  ■  Vi,  j  Vi.j  "yi-ljj 


m-53 


"0  '  '  >  o  • 

x.,.=x..+{x..-x.,.)  =  Zx..-x. 

i-lij  iiJ  i,J  i+lij  i.J  i+l.J 


ni-54 


With  these  initial  values,  the  iteration  can  be  ordered 
as  follows: 

From  111-49 


>-l 


if  f  -UK  f'y",  .wa  (x'V  +b 

J  1»J  1»J  J 


UI-55 


From  III-47 


h  II  h-1  h-1  ,,2 

(m  .  )  =  (x  .  .  )  -  (k.  .)  y 

i»J  hi  hi 


,  ^  II 


III-56 
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From  III-50 


II. h 


(K  .) 

i-l.j  ‘.J 


h-l 


h-l 


'i-l.jl"  )  ) 


h-l 
+  c 


III-37 


"2  ^ 


"h-l  1  1 

^-1  i  )  +  c.  ,  •) 

‘  1-1, jj  i.J 


Finally,  from  III-5Z 


2 


h 

)  + 


b. 


i-1.  j 


(x")^ 


+  c 


j  J 


m-58 


Thus,  with  the  initial  values  given  by  III-53  and  III-54,  equations 
III-55  through  III-58  are  iterated  until  convergence  is  obtained.  The 
criteria  for  this  will  be: 


where 

^  rn,  X  and  Q  is 


<  Q 


the  precision  required. 

The  application  of  the  above  technique  was  not  possible  in  this 
program,  due  to  limitation  of  time  and  funds.  It  is  felt,  however,  that 
the  analytical  procedures  discussed  in  this  section  are  valid,  and  that 
further  effort  in  the  application  of  these  techniqie  s  in  conjunction  with 
phenomenological  insight  gained  in  the  experimental  programs  would 
lead  to  the  accomplishment  of  the  initial  objectives  of  the  theoretical 
approach  to  the  parachute  problem. 
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SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY 


The  evolution  of  a  family  of  high  Mach  number,  self-inflating  para¬ 
chute  configurations  in  this  one  year  program  has  indicated  an  extension  of 
the  state  of  the  art  in  speed  capability  of  parachute  decelerators  by  a  factor 
of  about  three  over  that  capability  accepted  as  a  practical  limit  for  such 
decelerators  at  the  time  of  the  initiation  of  this  program.  This  family  pro¬ 
vides  also  a  drag  to  weight  ratio  advantage  of  a  factor  of  the  order  of  6  over 
any  other  known  deceleralor  capable  of  operation  in  this  Mach  number  range. 

Work  prior  to  this  program  had  led  to  the  development  of  parachutes 
capable  of  proper  performance  at  Mach  number  of  a  little  more  than  Z.  Re¬ 
quirements  for  applications  at  high  Mach  numbers  could  not  be  satisfied. 
Indeed,  it  was  not  even  well  established  that  it  was  feasible  to  develop  a 
parachute  for  applications  much  above  Mach  number  2.  The  maximum  Mach 
number  at  which  satisfactory  canopy  performances  were  exhibited  in  the  final 
(Phase  II)  wind  tunnel  test  program,  was  limited  only  by  the  Phase  II  wind 
tunnel  capability  (M  =  4.  65).  While  this  report  was  being  written,  a  test 
program  was  conducted  by  the  contracting  agency  in  which  configurations 
designed  during  this  program  were  successfully  tested  to  M  =  6.  This  also 
represented  the  speed  limitation  of  the  wind  tunnel  used  for  these  tests  (the 
Phase  I  test  facility.  Tunnel  A,  VKF,  AEDC), 

Due  to  the  successful  performance  of  the  high  speed  parachute  designs 
evolved  in  this  program  a  large  amount  of  data  was  gathered,  however,  since 
the  test  programs  (particularly  Phase  II)  were  the  first  to  be  conducted  at 
high  Mach  numbers  with  successfully  performing  configurations,  much  of  the 
data  serve  only  to  define  directions  of  effort  in  future  programs. 

A.  Conclusions 


The  following  specific  conclusions  were  drawn  during  the  conduct  of 
this  program; 

(1)  Although  several  basic  geometric  shapes  were  found  to  perform 
satisfactorily  in  the  high  Mach  number  environments,  those  configura¬ 
tions  which  consisted  of  a  cone  frustrum  shape,  of  high  inlet  diameter 
to  maximum  diameter  ratio  (0.  9)  and  with  a  10  degree  inlet  cone  angle, 
combined  with  a  flat  roof  appeared  to  be  superior.  In  these  designs, 
the  use  of  a  low  total  porosity  (8  to  18  percent),  a  low  porosity  (2  per¬ 
cent)  inlet  cone,  and  a  high  porosity  roof  (14  to  34  percent)  resulted  in 
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canopy  configurations  (C-6,  C-7,  CR-2)  which  performed  in  an  out¬ 
standing  manner  throughout  the  tested  speed  range.  Hereinafter  these 
configurations  employing  minor  variations  in  the  basic  design  charac¬ 
teristics  described  above  are  referred  to  as  the  Hyperflo  family  of 
high  performance  parachutes. 

(2)  The  use  of  a  relatively  low  total  canopy  porosity  in  the  Hyperflo 
designs  provides  high  over -all  inflation  characteristics,  in  addition  the 
very  low  porosity  of  the  canopy  inlet  cone  leads  to  high  stability  charac¬ 
teristics  for  these  configurations,  and  high  inflation  characteristics  in 
the  inlet  region.  The  cone  frustrum  shape  of  the  canopy  minimizes  the 
total  area,  and,  hence  the  total  weight  and  bulk  as  compared  for  example, 
to  a  canopy  which  ’.s  designed  with  a  preshaped  hemispherical  roof.  The 
latter  design  results  in  twice  the  roof  area  for  a  given  frontal  area.  The 
necessary  porosity  to  satisfy  mass  flow  and  stability  requirements  in 

the  Hyperflo  configuration  is  provided  by  the  high  porosity  roof. 

(3)  Observation  of  these  canopies  under  subsonic  conditions  during 
wind  tunnel  shutdown  has  indicated  excellent  inflation  and  stability 
characteristics.  It  is  therefore,  felt  that  successful  operation  of  the 
high  speed  canopies  will  be  assured  w'ith  perhaps  minor  variations  in 
the  design  parameters  for  subsonic  speed  applications. 

(4)  Although  relatively  low  total  porosities  were  used  in  the  design  of 
the  test  models,  as  compared  to  previously  applied  (conventional) 
parachute  designs,  wide  ranges  of  porosity  (8  to  18  percent)  were  used 
with  successful  results. 

(5)  The  canopy  location  aft  of  the  forebody  appeared  to  significantly 
affect  the  performance  of  all  configurations  in  the  high  speed  ranges. 

The  most  favorable  location  appeared  to  be  in  the  close-in,  or  low 

x/d  positions,  (6-8  forebody  diameters).  This  would  result  in  a  signi¬ 
ficant  reduction  in  the  over-all  parachute  system  weight,  since  the 
riser  length  will  be  minimized.  The  weight  of  a  parachute  riser  fre¬ 
quently  represents  a  large  percentage  of  the  total  parachute  system 
weight. 

(6)  The  effects  of  diameter  ratio  of  the  canopy  to  that  of  the  fore¬ 
body  was  not  clearly  established,  although  there  were  indications  of  a 
slight  reduction  in  canopy  performance  at  a  ratio  of  3.  33  compared  to 
a  ratio  of  2.  5. 

(7)  In  the  event  of  application  of  the  high  performance  parachutes 
to  an  unsymmetrical  forebody,  it  is  expected  that  satisfactory 
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parachute  per fonn.iiice  will  be  attained.  If  a  symmetrical  wake  is 
essential  to  satisfactory  canopy  performance,  the  location  of  the  canopy 
further  downstream  in  such  applications  should  be  satisfactory.  A 
pre\nous  application  of  a  coiiventional  parachute  decelerator  to  an  un- 
symmetrical  forebody  resulted  in  satisfactory  performance  at  M  =  2.2. 

(H)  The  effects  of  dpiamic  pressure  on  canopy  performance  are  not 
definitely  established,  .although  some  data  (»athered  in  this  program 
appear  to  indicate  high  drag  coefficients  witli  low  dynamic  pressure. 

The  effects  of  dynamic  pressure,  or  of  Reynolds  number  were  not 
investigated  sufficiently  in  this  program  to  establish  a  definite  correla¬ 
tion  with  canopy  performance. 

(9)  Several  roof  designs  were  employed  in  the  models  tested  in  this 
program.  No  significant  effect  on  canopy  performance  could  be  attri¬ 
buted  to  the  flat  roof  design  variation.  All  canopies  which  were  con¬ 
structed  using  the  basic  Hyperflo  design  geometry  performed  in  a 
comparable  manner  when  designed  with  equivalent  total  porosities. 

This  provides  a  wide  latitude  of  techniques  to  be  used  in  the  fabrication 
of  full  scale  canopies. 

(10)  Rigorous  analytical  procedures  by  means  of  which  canopy  per¬ 
formance  at  supersonic  speeds  might  be  predicted  was  completed  but 
not  applied  in  this  program.  The  validity  of  the  formulation  appeared 
to  be  established.  The  results  of  the  program,  however,  indicate  that 
relatively  non-rigorous  procedures  may  be  fruitful.  Careful  data 
analysis  coupled  with  relatively  non-rigorous  analytical  procedures  led 
to  the  indication  of  configurational  requirements  for  satisfactory  high 
Mach  number  parachute  performance  and  some  of  the  associated  air 
flow  phenomena.  Although  the  non-rigorous  procedures  applied  in  this 
program  resulted  in  successful  supersonic  canopy  performance,  a 
positive  prediction  of  supersonic  parachute  requirements  will  be  estab¬ 
lished  only  by  rigorous  analytical  methods. 

B.  Recommendations 


1.  Experimental  Programs 

a.  A  study  of  the  effects  of  dynamic  pressure  on  canopy  per¬ 
formance  in  the  high  Mach  number  range  is  essential,  due  to 
possible  dynamic  pressure  effects  observed  in  the  Phase  II  test 
program. 
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b.  The  effects  of  diameter  ratio  of  parachute  to  forebody  should 
be  investigated  to  determine  if  a  practical  limit  to  this  ratioexists. 
Although  two  such  ratios  were  investigated  in  the  final  test  pro¬ 
gram  (2.  5  and  3.  33),  larger  diameter  ratios  should  be  investiga¬ 
ted. 


c.  Further  investigation  to  determine  the  upper  and  lov/er 
limits  of  total  porosity,  within  which  successful  canopy  perform¬ 
ance  is  assured,  should  be  conducted. 

2.  Anal^^ical  Study 

Inherent  complexities  in  the  formulation,  and  necessary  mechani¬ 
zation  of  the  flow  field  solution  prevented  a  useful  contribution  of  this 
approach  to  the  canopy  designs  achieved  in  this  program.  It  is  felt  that 
the  feasibility  and  validity  of  the  analytical  procedures  discussed  in 
Section  III  appear  to  be  established,  however,  further  study  and  ex¬ 
perience  with  the  application  of  this  procedure,  by  matching  test  data 
and  incorporating  observed  phenomena,  should  result  in  an  invaluable 
contribution  to  future  parachute  canopy  designs,  and  the  further  under¬ 
standing  of  th'-se  phenomena. 
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APPENDIX  I 


THE  PH.ASE  1  TEST  FACILITY  DESCRIPTION 

The  VKF  Tunnel  A,  AEDC  Tennessee  is  a  continuous,  closed  circuit, 
variable  density  wind  tunnel  with  a  Mach  number  range  of  1.  5  to  6.  0.  The 
tunnel  is  served  by  the  main  compressor  system  which  provides  a  wide 
range  of  mass  flows  and  pressure  ratios  up  to  a  maximum  pressure  of  200 
PSIA. 


The  test  section  is  40  by  40  inches,  with  its  effective  test  chambers 
approximately  centered  at  the  downstream  end  of  the  flexible  nozzle  plates. 
Three-pane  windows  are  located  in  the  test  section  side  walls. 

A  variable -geometry,  five-hinge  diffuser  is  located  downstream  of  the 
test  sc'tion. 

An  ERA-1102  digital  computer  is  shared  by  all  VKF  tunnels,  with 
individual  tunnels  equipped  with  a  digital  scanner,  raw-data  punch,  and 
related  output  equipment. 

A  double-pass  Schlieren  system  with  a  35  inch  diameter  optical  path 
is  provided  for  observation  of  flow  fields  through  any  of  the  three  test 
section  windows. 
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APPENDIX  II 


THE  PHASE  I,  AND  la  TEST  FACILITY  DESCRIPTION 


The  Cook  Technological  Center  wind  tunnel  Is  a  continuous  flow  open 
circuit  atmosphere  induction  type  with  either  fixed  or  variable  flow  nozzles. 
The  test  sections  are  of  the  two-dimensional  type  with  airflow  available  by 
adjustment  of  two  General  Electric  J--J7  turbojet  engines.  The  facility  was 
designed  to  provide  continuous  variation  of  Mach  numbers  from  low  super¬ 
sonic  to  a  Mach  number  of  approximately  2,  7  in  a  test  section  22  inches  by 
24  inches. 

Flow  visualization  is  accomplished  by  the  utilization  of  a  Schlieren 
optical  system  which  is  adjustable  by  sliding  on  rails  to  permit  rv  cording 
access  to  all  sections  of  the  test  chambers.  Two  window  ports  in  the  side 
w'alls  are  utilized  to  photograph  shock  wave  patterns  at  different  sections 
of  the  te  St  area. 
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APPENDIX  III 


4 


THE  PHASE  II  TEST  FACILITY  DESCRIPTION 

I 

The  Langley  Large  Supersonic  Tunnel,  Langley  Field,  Virginia  is 
capable  of  test  Mach  numbers  from  2.  3  to  4.65,  with  a  maximum  stagnation 
pressure  of  150  PSIA.  The  test  section  of  this  facility  is  4  feet  high,  4  feet 
wide,  and  7  feet  long,  and  will  permit  variat  ion  of  Mach  number  at  any 
desired  increment  throughout  its  range  with  tunnel  operating. 

Both  stagnation  pressure  and  stagnation  temperature  may  be  controlled 
independe  ntly. 

Tlie  tunnel  window  has  a  field  of  59  inches  by  48  inches,  and  is  made 
up  of  nine  strips  of  optical  plate  glass,  each  5-1/2  inches  wide,  and  separa¬ 
ted  from  each  other  by  1-1/4  inches  of  supporting  structure. 

The  Schlieren  system  has  a  49  inch  diameter  field,  and  adjustable  so 
as  to  allow  complete  coverage  of  the  59  inch  window. 

A  semi-automatic  force  data  readout  system  provides  tabulated  raw 
data  and  IBM  punch  card  storage  of  raw  data  concurrent  with  tunnel  opera¬ 
tion. 

t 
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